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ABSTRACT 

Aims. We present an analysis of the relation between the star formation rate (SFR) surface density (Esfr) and mass surface density 
of molecular gas (Sh 2 ), commonly referred to as the Kennicutt-Schmidt (K-S) relation, on its intrinsic spatial scale, i.e. the size of 
giant molecular clouds (—10—150 pc), in the central, high-density regions of four nearby low-luminosity active galactic nuclei (AGN). 
These are AGN extracted from the NUclei of GAlaxies (NUGA) survey. This study investigates the correlations and slopes of the 
K-S relation, as a function of spatial resolution and of the dilferent 12 CO emission lines used to trace Sh 2 , and tests its validity in the 
high-density central regions of spiral galaxies. 

Methods. We used interferometric IRAM 12 CO(1-0) and l2 CO(2-l) and SMA 12 CO(3-2) emission line maps to derive and 
HS T-Wa images to estimate £ S fr- 

Results. Each galaxy is characterized by a distinct molecular SF relation on spatial scales between 20 to 200 pc. The K-S relations can 
be sublinear, but also superlinear, with slopes ranging from ~0.5 to ~ 1.3; slopes are generally superlinear on spatial scales > 100 pc and 
sublinear on smaller scales. Depletion times range from ~1 and 2 Gyr, which is compatible with results for nearby normal galaxies. 
These findings are valid independently of which transition - l2 CO(l-0), l2 CO(2-l), or 12 CO(3-2) - is used to derive . Because 
of either star-formation feedback, the lifetime of clouds, turbulent cascade, or magnetic fields, the K-S relation might be expected to 
degrade on small spatial scales (<100 pc). However, we find no clear evidence of this, even on scales as small as ~20 pc, and this 
might be because of the higher density of GMCs in galaxy centers that have to resist higher shear forces. The proportionality between 
E h , and E SF r found between 10 and 100 M Q pc~ 2 is valid even at high densities, ~10 3 M 0 pc -2 . However, by adopting a common 
CO-to-H 2 conversion factor (a C o), the central regions of the NUGA galaxies have higher E S fr for a given gas column than those 
expected from the models, with a behavior that lies between the mergers or high-redshift starburst systems and the more quiescent 
star-forming galaxies, assuming that the first ones require a lower value of a co- 
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1. Introduction 


The relationship between gas and star formation (SF) in galaxies 
plays a key role in galaxy evolution. It constrains how efficiently 
galaxies turn their gas into stars an d also serves as essential input 
to simulations and models (e.g., Bpissjer^taL 2003UTanetalJ 


1999fc ISnringel & Hernauist 1200? Krumholz & McKeel 120051) . 

Despite the importance of this relationship, the processes respon¬ 
sible for the conversion of gas into stars in various galactic envi¬ 
ronments are still poorly understood. 

More than fifty years ago, ISchmidtl dl959t 119631) suggested 
that the star formation rate (SFR) and the gas content in galaxies 
are related by 


2sfr=A^ 


( 1 ) 


where Esfr and E gas are the SFR surface density and the gas sur¬ 
face density, respectively, A is the normalization constant repre¬ 
senting the efficiency of the processes regulating gas-stars con¬ 
version, and Nfa the power relation index. The gas can be atomic 
(Hi), molecular (H 2 ), or a combination of both (H 1 +H 2 ). 

Although the atomic phase of the interstellar medium (ISM) 
is directly traced by the Hi emission line at 21cm, indirect 
approaches are needed to estimate the distribution of H 2 . The 
molecular hydrogen indeed lacks a dipole moment and typical 


temperatures (15-25 K), in giant molecular clouds (GMCs) are 
too low to excite quadrupole or vibrational transitions. For these 
reasons, the carbon monoxide (CO) emission lines are the most 
straightforward and reliable tracer of H 2 in galaxies. CO is rela¬ 
tively bright and its ability to trace the bulk distribution of H 2 has 
been confirmed by comp arisons with gamma ravs (iLebrun et alJ 


^1983^Shxm geLglJ ^988) and dust emission ( Desert et alJ _ 
Dam e et alJ 20011) . Comparing the atomic gas density with the 


number of representative young stellar objects in the solar neigh¬ 
borhood, [Schmjdtj ( 1959 ) derived a power relation index iV»2 
in equation (©. Values of N ranging appr oximately from 1.5 to 2 
were further confirmed bv lGuibert et al. ( 1978 ). using more pre¬ 
cise data on the radial and vertical distributions of the interstellar 
gas and a variety of young objects as tracers of SFR. 


Later. iKennicutd (il998alb!) found that for disk-averaged sur¬ 
face densities, both normal star-forming and starburst galaxies 
follow equation © with a power relation index of N ^ 1.4 for 
total (H 1 +H 2 ) gas. This correlation between gas and SFR surface 
density is commonly referred to as the Kennicutt-Schmidt rela¬ 
tion (hereafter the K-S relation). Such a relation is in principle 
consistent with large scale gravitati onal instability being the ma - 
jor driver of dense cloud formation (IOuirklll972llMadoral 19771) . 
Although spatially unresolved studies of Hi, CO, and SFR are 
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useful for characterizing global disk properties, understanding 
the mechanisms behind the SF requires resolved measurements. 

It is now possible to study the K-S relation at sub-kpc scales, 
approaching the intrinsic scale of SF, i.e. the size of GMCs 
(~ 10—150 pc, e.g.. ISolomon et all 1 19871) . Several factors have 
contributed, including the explosion in multiwavelength data for 
nearby galaxies: for example, UV-GALEX: Gil de Paz et al. 
2007; infrared -Spitzer and H a-HST: SINGS and LVL surveys, 
Kennicutt et al. (2003), Dale et al. (2009); far-infrared-Herschel: 
KINGFISH, Kennicutt et al. (2011); 12 CO: BIMA SONG, Heifer 
et al. (2003); HERACLES-IRAM 30m, Leroy et al. (2009); Hi: 
VLA THINGS survey, Walter et al. (2008). Also technical im¬ 
provements, especially the construction of millimeter interfer¬ 
ometers, now allow higher resolution imaging of 12 CO in galax¬ 
ies. 


There is now substantial evidence that the molecular gas 
is well-correlated with the SFR tracers over several orders of 
magnitude, but mostly for regions where H 2 makes up the 


majority of the neutral gas, Xh 1 L Xhi (e.g., Wong & Blitz 

2002 

; Komugi et al. 

2005; Kennicutt et al. 2007; Thilkeretal. 

2007 

Schuster et al. 

2007; Bigiel et al.l 200i§ Lerov et al. 

2008; 

Blanc et al. 2009; 

Onoderaetal. 2010; Verlev et al. 

2010; 

Rahman et al. 2011 

Momose et al. 2013». This is because 


GMCs, the major reservoirs of molecular gas, are the sites of 
most SF in our Galaxy and other galaxies. Their properties in¬ 
deed set the initial conditions for protostellar collapse and may 
play a role in determinin g the stellar initial mass function (IMF, 
iMcKee & Qstrikeil2007l ). Moreover, it has long been known that 
the spatial distributio n of 12 CO emission follows that of the stel- 


lar light and Her (e.g.. Young & Scoville 1982; 

Scoville & Young 

1983; Solomon et al. 1983; Tacconi & Young 

1990). The lack of 


a clear correlation between Shi and X SFR inside galaxy disks 
(e.g.. iBigiel et alJl2008l) offers circumstantial evidence that SF 
remains coupled to the molecular, rather than total gas surface 


density (Shi ihi) e ven where Hi makes up most of the ISM. 
[Bigiel et all ( 2008 ) find that Shi saturates at surface density of 
~9 M 0 pc 2 , and gas in excess of this value is in the molecular 
phase of spirals and Hi-dominated galaxies. 

A crucial parameter in the study of the resolved K-S relation 
is the choice of the tracer for Hi and SFR. Most of the studies 
for deriving the H 2 surface density to study the molecular SF 
relation are based on the 12 CO(1-0) (e.g., iRahman et all [20171 
ISaintonge et alJ 2011J;|MQmpseW_ aj]|2QI3|) or 12 CO( 2-l) emis¬ 
sion lines (e.g., Bigiel et all 20081 : iLerov et~ak 20081), and more 
rarely on the 12 CO(3-2) transition (e.g., iKomugi et al.1120071 : 
Wils on et al.ll2012l) . Also the dependence of the SFR on dense 
molecular gas mass has been explored with molecules such as 
HCN and HCO + (e g.. iGao & Solomonl 12004 IWuetaD 120051: 
iGarcfa-Burillo et alJl2012h. 

Over the past thirty years, extensive efforts have been made 
to der i ve plausible SFRs f or external galaxies (see iKennicutl 
Il998at ICalzetti et afl 1201 (I for reviews). Optical SFR tracers, 
such as Her, often suffer from dust extinction, which can change 
dramatically from location to location. ICalzetti et alJ ( 2007 ) find 
Ay ~ 2.2 mag in typical extragalactic Hu regions and dense 
star-forming regions can be completely obscured with Ay that 
can reach ~6 mag ('e.g JScoville et alJl2001h . In addition, the Ay 
value is not always known for a given galaxy. However, IR space 
facilities (Spitzer and Herschel) and UV (GALEX) have made it 
possible to also image the details of dust-obscured SF. 

Different SFR tracers probe different time scales and hence 
the SF history of any particular galaxy. Her emission traces ion¬ 
ized gas by massive (M > 10 M e ) stars on a timescale of 
< 20 Myr. The far-UV (FUV) luminosity corresponds to rela¬ 


tively older (< 100 Myr), less massive (> 5 M e ) stellar popu¬ 
lations. The MIR 24 pm emission mostly traces reprocessed ra¬ 
diation of newborn (few Myr) OB stellar associations embed¬ 
ded within the parental molecular clouds. Although star clusters 
emerge from their natal clouds in less than 1 Myr, they remain 
associated with it on a much longer time scale, ~ 10-30 Myr, 
which is the ti me scale assoc i ated w ith t he 24 pm emis s ion as 
a SFR tracer. iKennicutt et al.l (1 20071) and ICalzetti et al.l {2007) 
have justified the feasibility of correcting the number of ionizing 
photons, as traced by the Hu recombination line, for the effects 
of the dust extinction by adding a weighted component from 
the Spit zer MIPS 24 pm lu mino sity in indiv i dual s tar-forming 
regions. ILerov et all (120081) and IBigiel etafl d2008l) have pro¬ 
posed another composite SFR tracer that corrects dust attenu¬ 
ation of far-UV surface brightness using, also in this case, the 
Spitzer 24pm emission but with different weights. These com¬ 
posite SFR tracers, Ha+24pm and FUV+24/rm, have been ex¬ 
tensively used for a large number of nearby galaxies (see refer¬ 
ences above). 


Spatially resolved K-S relation st u dies p ublished in 


the 

past decade (e.g., Kennicutt et al. 2007| Bigie 

et al. 

2008; Lerov et al. 2008; Blanc et al. 2009; Bigiel et al. 

2010; 

Eales et al. 2010: Verlev et al. 2010; Bigiel et al. 2011; Liu et al. 

2011 

; Rahman et al 201 , 2012; Schruba et al. 2011; Ford et al. 

2013 

; Viaene et al. 2014) found a wide spread in the value of 


the power relation index (N ~ 0.6 - 3) both within and among 
galaxies. For a comprehensive review of the most recent SF re¬ 
lation studies on sub-kpc scales, we refer to IKennicutt & Evansl 
(120121) . The wide range of the value of N may be intrinsic, sug¬ 
gesting that different SF ‘laws’ exist, and may contain valuable 
astrophysical information. Alternatively, it may be partially or 
entirely due to the adop ted physical scale (e.g., ISchruba et al.l 
2(H0l;ICalzettLetalJl2Q1 2|), the choice of the molecular gas (e.g., 
Krumholz & Thompsonll2007t iNaravanan et al] F2008) and SFR 
tracer, the type of galaxy under investigation, data sampling, 
and th e fitting method used (e.g., Blan c et akll2009l : IShettv et all 
120131) . In particular, different SFR tracers and spatial scales ef¬ 
fectively sample different timescales, so that a galaxy’s SF his¬ 
tory can play a role in determining the results of the measure¬ 
ments. It is also possible that these differences correspond to 
a spectrum of physical mechanisms present in a wide range 
of environments. High shear in galactic bars, harassment in a 
dense galaxy cluster, and galaxy mergers and interactions have 
the potential to eithe r dampen or enhance the SF process (e.g., 
ICasasola et alJ 12004 IZhou et all 1201 lb lLanz et all 12013l) . All 
these environmental processes at work, on both galactic and ex¬ 
tragalactic scales, suggest therefore that there is no universal SF 
relation in the Local Universe. 

In addition to the slope of the empirical power relation re¬ 
lation between X gas and Xsfr, the other crucial parameter in SF 
studies is the molecular gas depletion time, defined as the time 
needed for the present SFR to consume the existing molecular 
gas reservoir, Td ep i = £ gas /EsFR, i- e - the inverse of the star for¬ 
mation efficiency (SFE). One interpretation of the linear slope 
for the K-S relation is an approximately constant Td ep i, with an 
average Td ep i of about 2 Gyr in normal spirals (e.g.. ILerov et al.l 
120131) . As for the power relation index, several factors compli¬ 
cate the interpretation of Td ep i- Moreover, while a linear slope 
(and constant Td ep i) describes the global average scaling rela¬ 
tion in the star-forming gala xies well, individual galaxies devi¬ 
ate fr om a single Td ep i (e.g., Saintong e et al.ll20Tlt ILerov et all 
120131) . 

The objective of this paper is to explore the molecular SF re¬ 
lation in the inner (~1 kpc), high-density regions of four nearby 
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(D < 20 Mpc), low-luminosity active galactic nuclei (LLAGN) 
on the spatial scale of ~20-200 pc through a pixel-by-pixel anal¬ 
ysis of the available maps. We do not discuss the correlation of 
SFR with the atomic component of the gas, since the gas phase is 
predominantly molecular in the central regions we study. These 
galaxies were originally part of the NUclei of GAlaxies (NUGA) 
survey carried out at the Plateau de Bure Interferometer (PdBI, 
iGarcfa-Burillo et aI1l2003l) . The good spatial resolution offered 
by 12 CO NUGA observations enables us to probe GMC spatial 
scales around LLAGN, and better understand the SF process. 
Here, we investigate the correlations and slopes of the SF rela¬ 
tion as a function of spatial resolution and of various 12 CO tran¬ 
sitions (1-0, 2-1, 3-2) for deriving Zh 2 , and the validity of the 
K-S relation in regions with high molecular gas surface densi¬ 
ties. 

The NUGA sample is ideal for such a study for several rea¬ 
sons: The proximity of the galaxies combined with the good spa¬ 
tial resolution and sensitivity afforded by PdBI and SMA and the 
high-resolution SFR tracers from the Hubble Space Telescope 
( HST ) let us probe the K-S relation on fine spatial scales that up 
to now have only been examined in Local Group galaxies. More¬ 
over, the physical conditions in the NUGA LLAGN are more 
extreme than those found in typical spiral disks. We can thus as¬ 
sess how well the warm dense gas at high column densities in 
AGN circumnuclear regions can form stars. These findings set 
the stage for what will be possible with ALMA data. 

This paper is organized as follows. We describe the sample 
selection in Section [2] and the data we used, their treatment, the 
procedure for deriving Eh, and Xsfr maps from the original im¬ 
ages, and the applied fitting method in Sect. [3 In Sects. [4] and 0 
we show the results on the observed relationships between Zh 2 
and Xsfr for individual galaxies and the whole sample, respec¬ 
tively. In Sect. [6] we stress the caveats and uncertainties associ¬ 
ated with the present study, and finally in Sect. [7] we summarize 
our findings and give our conclusions. 

2. The sample 

The sample presented in this paper consists of four LLAGN se¬ 
lected from the NUGA survey. The NUGA project is an IRAM 
Plateau de Bure Interferometer (PdBI) and 30 m single-dish sur¬ 
vey of nearby LLAGN with the aims (i) of mapping, at the an¬ 
gular resolution of ~0.5-2" and sensitivity of ~2-4 mJybeanR 1 , 
the distribution and dynamics of the molecular gas through the 
two lowest emission lines of the 12 CO in the inner kpc of the 
galaxies, and (ii) of studying the different mechanisms for gas 
fueling of LLAGN. Each galaxy of the core NUGA sample (12) 
has been studied on a case-by-case basis. 

In this paper, we present the results of the study of the molec¬ 
ular gas spatially-resolved SF relation in the following NUGA 
galaxies: NGC 3627, NGC 4569, NGC 4579, and NGC 4826. 
This NUGA subsample offers heterogeneity in terms of nuclear 
activity, distance, detection or not of gas inflow, gas morphol¬ 
ogy, and surrounding environment. Since the sample only con¬ 
sists of four galaxies, the wide range of galaxy properties cannot 
be used to infer statistical conclusions. It is worthwhile, how¬ 
ever, looking at the different conditions where the SF relation is 
studied in the selected NUGA subsample. These properties are 
collected in Table |T] In this table. Col. (1) indicates the galaxy 
name. Cols. (2) and (3) the coordinates of the galaxy dynami¬ 
cal center derived from NUGA IRAM 12 CO observations (see 
later), Col. (4) the morphological t ype from the Third Referenc e 
Catalog of Bright Galaxies (RC3. Ide Vaucouleurs et alJ 1 199 lh . 
Col. (5) the nuclear activity, Col. (6) the distance ( D ), Col. (7) 


the inclination (/), Col. (8) the position angle (PA), Col. (9) the 
prevalent molecular gas morphology as detected from NUGA 
observations. Col. (10) the identified molecular gas motion (i.e., 
inflow). Col. (11) the surrounding environment, and Col. (12) 
the NUGA references. 


3. The data, derived parameters, and fitting method 

This study is based on measurements of the molecular gas and 
SFR surface densities, Xh 2 and Xsfr, and relies on the existence 
of multiwavelength datasets. We used l2 CO(l-0), 12 CO(2-l), 
and i 2 CO(3-2) line intensity maps to derive the surface densities 
of the molecular gas, and HS T Her (6563 A) emission images to 
estimate the surface densities of SF. In this section we describe 
these datasets. 

3.1. IRAM l2 CO(1-0) and n CO(2-1) observations 

The 12 CO(1-0) and l2 CO(2-l) line intensity maps are part of the 
NUGA survey. NUGA observations have been carried out with 
six antennas of the PdBI in the ABCD configuration of the ar¬ 
ray and with the IRAM 30m single-dish telescope. 12 CO images 
were reconstructed using the s tandard IRAM/GILDAffl soft¬ 
ware dGui lotea u & Lucas 2000), following the prescriptions de¬ 
scribed in Garc fa-Burill o et a7Td2003l) . and restored with Gaus¬ 
sian beams. The beam sizes are typically %2" for 12 CO(1-0) 
and SU" for 12 CO(2-l). We used natural and uniform weight¬ 
ings to generate 12 CO(1-0) and 12 CO(2-l) maps, respectively. 
This allowed us to maximize the flux recovered in 12 CO(1-0) 
and optimize the spatial resolution in 12 CO(2-l). For NGC 3627 
and NGC 4579, 30m observations were used to compute the 
short spacings and complete the interferometric measurements, 
whereas for NGC 4569 and NGC 4826, we used 12 CO maps ob¬ 
tained only from PdBI observations. 

For each NUGA galaxy we derived the position of the AGN 
by assuming that this coincides with the dynamical center of the 
galaxy (see Table[l}. This choice maximizes the symmetry of the 
global velocity field derived from the two 12 CO transitions. De¬ 
tails on the determination of the AGN position/dynamical cen¬ 
ter are discussed on a case-by-case NUGA papers. Properties of 
IRAM I2 CO maps are collected in Table [2] and described in de¬ 
tail in Sect. 13.51 

3.2. SMA l2 CO(3-2) observations 

For NGC 4569 and NGC 4826, we also have 12 CO(3-2) line in¬ 
tensity maps obtained with the Submillimeter Array (SMA) in its 
compact configuration with seven working antennas. 12 CO(3-2) 
single-dish observations carried out with the 15 m J ames Clerk 
Maxw ell Telescope (JC MT) and pu blished by IWilson et alJ 
(120091) for NGC 4569 and llsraell (l2009lt for NGC 4826 were used 
to compute and add the missing short spacings to the SMA data. 
The beam sizes of 12 CO(3-2) maps are comparable to those of 
the IRAM 12 CO(1-0) maps, i.e. ~2" (see Table O. Details on 
these observations are contained in lBoone et al.1 (1201 ll) . 

3.3. HST Her emission-line images 

We retrieved Her emission-line images from the Hubble Legacy 
Archive (HLA) that makes HST WFPC2 observations of our 
galaxy sample available. In these images, the Her emission line 

1 http://www.iram.fr/IRAMFR/GILDAS/ 
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Table 1: Properties of the galaxy sample. 


Galaxy 

042000 

dj2000 

RC3 type 

Nuclear act. 

D 

i 

PA 

Gas morph. 

Gas mot. 

Envir. 

NUGA Ref. Ih) 


|"h m sj 

[° ' "1 



[Mpc] 

[deg] 

[deg] 





(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(ID 

(12) 

NGC 3627 

11 20 15.02 

12 59 29.5 

SAB(s)b 

LINER/Seyfert 2 

10.2 (1"=49 pc) 

61 

178 

Bar 

Inflow 

Interacting 

(1) 

NGC 4569 

12 36 49.80 

13 09 46.3 

SAB(rs)ab 

Transition 2 

16.8 (1"=81 pc) 

62 

15 

Bar 

Inflow 

Isolated 

(2,3) 

NGC 4579 

12 37 43.52 

11 49 05.5 

SAB(rs)b 

LINER/S y 1.9 

19.8 (l"=96pc) 

36 

95 

Spiral arms 

Inflow 

Isolated 

(4,5) 

NGC 4826 

12 56 43.64 

21 40 59.3 

(R)SA(rs)ab 

LINER 

5.4 (1"=26 pc) 

54 

112 

Disk 

- 

Isolated 

(3,4,6) 


No tes. (a) Molecular g as m orphology as reveal ed fr om NUGAIRAM 12 CO(l-Q ) ob servations. For 12 CO ( 2-l) a nd l 2 CO(3-2) morpholog i es we refer to NUGA papers. (i>) (ll lCasasola et all d201 ll) : 
(2) lBoone et all d2007t) : Gl lBoone et all d201 lh : (4) lGarci'a-Biirillo et alJd2005h : (5l lGarcfa-Burillo et all d2009h : (bl lGarci'a-Burillo et alld2003l) . 
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Table 2: Basic information on 12 CO NUGA dataset. 


Galaxy 

IRAM instrument 

12 CO line 

Beam 

lcr 

f 10, F 21, F 32 

Rn. R 32 




[" x ", °] 

[Jy beam" 1 km s" 1 ] 

[K(Jy beam" 1 )" 1 ] 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

NGC 3627 

PdBI+30m 

d-0) 

2.05 x 1.26, 30 

0.16 

35.55 

- 


PdBI+30m 

(2-1) 

0.86x0.58,23 

0.30 

46.32 

0.60 

NGC 4569 

PdBI 

(1-0) 

2.33 x 1.46, 27 

0.25 

27.00 

- 


PdBI 

(2-1) 

1.17x0.65, 198 

0.35 

30.24 

0.63 


SMA+JCMT 

(3-2) 

2.60 x 2.06, -36 

6.19 

1.92 

0.23 

NGC 4579 

PdBI+30m 

(1-0) 

2.03 x 1.28,26.03 

0.16 

35.34 

- 


PdBI+30m 

(2-1) 

0.98 x 0.60, -155.48 

0.32 

38.94 

0.90 

NGC 4826 

PdBI 

(1-0) 

2.53 x 1.80 

0.25 

20.17 

- 


PdBI 

(2-1) 

0.77 x 0.55 

0.36 

54.21 

0.67 


SMA+JCMT 

(3-2) 

2.58 x 1.94 

4.44 

2.05 

0.38 
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was observed through the narrow-band filters F656N or F658N, 
and the underlying continuum through F547N, F555W, and/or 
F814W (equivalent to narrow-band V, V, and /, respectively). 
The HS T images have a pixel size of 0'.'05 and a spatial resolu¬ 
tion of O'/l. 


The available maps are emission-line-only F[q-+[Nii] images 
containing emission from both Her at 6563 A and [Nii]/Tl 6548, 
6584 A. We removed the [Nn] contamination within the filter 
bandpass using average [N n]/Ha values available in literature. 
Then, we corrected the Her maps both for Galactic foreground 
extinction and internal extinction. We assumed the values of the 
/7-band Galactic foreground extinction ,4#(Gal) available from 
the literature (see Tabled and used the interstellar extinction 
curve by ICardelli et al.1 (1989). For three galaxies of our sam¬ 
ple (NGC 3627, NGC 4569, and NGC 4826), interna l extinction 
corrections were extracted from lCalzetti et all d2007l) using Pact 
images as a yardstick for calibrating the MIR emission. The hy¬ 
drogen emission lines trace the number of ionizing photons, and 
the Pact line (at 1.8756 pm) is only modestly affected by dust 
extinction. Because of its relative insensitivity to dust extinction 
(less than a factor 2 of correction for typica l extinction in nearby 
galaxies. Ay < 5 mag. ICalzetti et al. 12007), Paa is a nearly un¬ 
biased tracer of the current SFR ( Kennicuttll 1 998bl) . Among the 
various SFR calibrations, the linear combination of the observed 
Ha and the 24 pm emission is the one most tightly correlated 
with the extinction-corrected Paa emission (ICalzetti et al.l2007h . 
The most straightforward interpretation of this trend is that the 
24 pm emission traces the dust-obscured SF, while the observed 
Ha emission traces the unobscured one (Kennicutt et al.112007 ). 
Thus, the combination of the two recovers all the SF in a given 
region. This interpretation is confirmed by models (e.g., STAR- 
BURST99, 2005 update, Leitherer et al. 2005; Draine & Li 
2007), which also suggest that the trend is relatively independent 
of the characteristics of the underlying star-forming population. 
This implies that Ha/Paa ratio and the internal extinction de¬ 
rived from it, for instance Ay(int) in the V -baud as performed by 
ICalzetti et al.l (120071) . are good tracers of the internal extinction 
of a galaxy. 


For NGC 4579, rather than the Ha/Paa ratio, we used 
the internal extinction derived from the Balmer decrement by 
iHol ( 1999 ). One of the most reliable techniques for estimat¬ 
ing interstellar extinction is indeed to measure the flux ratio 
of two nebular Balmer emission lines such as Ha/H/j (i.e., the 
Balmer decrement). The determination of dust extinction from 
the Balmer decrement has been shown to be a very successful 
technique in the Local Universe since the first statistical work 
bv lKennicutt ( 1992 ). These results have been improved upon by 
the large amount of optical spectra provided by the Sloan Dig¬ 
ital Skj^SurveyA^DS^which were analyzed in t his context 
bv lBrinchmann etakl (1 2004) and lGarn & Best! (l2010h . The main 
limitation to applying the internal extinction correction derived 
from Balmer decrement is to neglect a possible completely ob¬ 
scured SFR component, hence to underestimate the total SFR. 
In any case, NGC 4579 has a Ay(int) value similar to those of 
the other galaxies, as shown in Table [2 This table collects the 
properties of the original Ha+[Nn] images and the values of the 
parameters used to obtain final Ha maps. In Tabled Col. (1) 
indicates the galaxy name. Col. (2) the HS T instrument. Col. 
(3) the ler noise of the background subtracted Ha+[Nn] images. 
Col. (4) the adopted [Nn]/Ha ratio. Col. (5) the fi-band Galactic 
foreground extinction |4/ ; (Gal)], and Col. (6) the V-band inter¬ 
nal extinction [Ay(int)]. 


3.4. Image treatment 

With the images described above, we constructed maps of In, 
and Zsfr to perform a pixel-by-pixel analysis of the spatially re¬ 
solved K-S relation. The procedures for obtaining the final maps 
are described later in Sections 13.51 and 13.61 Since we compared 
images with different properties and wide spreads in resolution, 
the first task was to convert them to a common alignment and 
resolution. 

All the original images ( 12 CO(1-0), 12 CO(2-l), l2 CO(3-2), 
and Ha) were centered on the dynamical centers of the galax¬ 
ies derived fromNUGA IRAM 12 CO observations (see TableQ}. 
The Ha maps were convolved to the resolution of the 12 CO (1- 
0, 2-1, and 3-2) maps with a Gaussian beam (on the sky); i.e., 
we do not account for the inclination of the galaxy. Then, all 
maps have been resampled to a pi xel size equal to the adopted 
12 CO resolution (see, for instance, Vutisalchavakul et al.ll2014l 
for a similar treatment of images). Since the spatial resolution 
and pixel size are equivalent, the pixels can be considered as 
roughly statistically independent, and there should be little cor¬ 
relation among them. Each pixel was thus treated as a single data 
point. Finally, we also convolved all maps to a common 200 pc 
resolution to be able to distinguish between the effects of differ¬ 
ent J transitions and spatial resolution. These procedures were 
performed by using IRAM/GILDAS and IRAffj softwares. 

We thus probed the gas and SFR surface densities in indi¬ 
vidual galaxies on physical scales ranging from 17 pc to 190 pc, 
as well as with a common resolution of 200pc. Such scales 
are smaller than those previously scrutinized ou t side the Local 
Group fe.g JBigiel et al.ll2008l; iLerov et al.ll2008L iRahman et al.l 
1201 lh and are superior in this respect to other kinds of analy¬ 
ses, s u ch as azimuthally averaged radial profile s (e.g JKennicuttl 
1989; iBigiel etafl 120081 : IRahman et al.1 l20flh or the aperture 
analysis encompassi ng star-forming regions and centering on H a 
emission peaks (e.g jKennicutt et aLll2007HBlanc et alJl2009l) . 


3.5. Molecular gas surface density maps 

The availability of 12 CO(1-0), 12 CO(2-l), and 12 CO(3-2) data 
offers In, maps at different resolutions and the possibility to 
compare the three lowest I2 CO transitions in terms of the K- 
S relation. We derived In, maps from 12 CO(1-0), l2 CO(2-l), 
and l2 CO(3-2) integrated intensity maps (Sco(i-O)- Scoa-i), 
S co(3-2)) by adopting a constant value for the Wo conversion 
factor , Xco = 2.2 x 10 20 cm 2 (K km s -1 ) -1 (ISolomon & BarretS 
19911) that corresponds.to aco = 3.5 M C ,, p c ~ (Kkm s” 1 )” 1 (e.g., 
Narayanan et al.l2012llBolatto et al.l2013l) . For 12 CO(1-0) emis¬ 
sion, the conversion to Ih, is 


Ih 2 - 3.5Sco(i-0)Eiocos/ 


( 2 ) 


where Ih, is in units of M Q pc~ 2 , Sco(i-O) i n JybeanD 1 kms -1 , 
Fio is the conversion factor from flux density to brightness tem¬ 
perature for the 12 CO(1-0) line in KfJybeam 1 ) “ 1 , and i is the 
galaxy inclination. For 12 CO(2-l) emission, the conversion to 
Ih, is derived from equation Q by defining R 2 l as the 12 CO(2- 
1)/ 12 CO( 1 —0) line ratio in temperature units: 


£ h . 


SCQ(2-1) 

R 2l 


F 2 i cos i 


(3) 


2 IRAF is the Image Reduction and Analysis Facility. IRAF is writ¬ 
ten and supported by the National Optical Astronomy Observatories 
(NOAO) in Tucson, Arizona. NOAO is operated by the Association of 
Universities for Research in Astronomy (AURA), Inc. under coopera¬ 
tive agreement with the National Science Foundation 
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Table 3: Basic information on the HS T -H a dataset. 


Galaxy 

HS T instrument 

lcr 

[Nii]/Ha- 

A b ( Gal) 

Ay(int) 



[erg s -1 cm -2 ] 


[mag] 

[mag] 

( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 

NGC 3627 

WFPC2(F658N) 

5.43 x 10 -iy 

0.54 (a) 

0.08 <a) 

1.94 <b) 

NGC 4569 

WFPC2(F656N) 

1.90 x 10 -18 

0.50 (b) 

0 . 20 (c) 

1.05 (b) 

NGC 4579 

WFPC2(F658N) 

9.49 x 10 -18 

0.62 (d) 

0.18 (c) 

2 . 10 (e) 

NGC 4826 

WFPC2(F656N) 

6.19 x 10 " 18 

0.51 (a) 

0.16 (a) 

1.97 <b) 


N otes. '"'iKennicutt et all {2008). These v alues are an average of the values fro m the maps of Burstein & Heiles ( 1982) and lSchlegel eTafl (TgH). 
m ICalzettietal]d2007h . ^ From NED. (J) lKennicutt etaD <2009h . w Ih 3 d!999h . 


where F 21 is the conversion factor from flux density to bright¬ 
ness temperature for the 12 CO(2-l) line. Consistently with equa¬ 
tion (0, for i 2 CO(3-2) emission the conversion to Xh 2 is derived 
by defining R 32 as the 12 CO(3-2)/ 12 CO(l-0) line ratio: 

= 3.5—-- F 32 cost (4) 

«32 

where F 32 is the conversion factor from flux density to bright¬ 
ness temperature for the 12 CO(3-2) line. Both Ri 1 and R 32 used 
in the equations {3} and ©, respectively, have been measured 
for each NUGA galaxy; we used the mean value of these mea¬ 
sured ratios to convert all transitions to CO(l-O). Equations (0. 
(0), and (0) define hydrogen surface densities; i.e., they do not 
include any contribution from helium. To scale our quoted sur¬ 
face densities to account for helium, they should be multiplied 
by a factor —1.36. 

The values needed to apply equations {2]), <0. and ({4} are 
listed in Tables 0 and 0 In Table 0 Col. (1) gives the galaxy 
name. Col. (2) the instruments used for observations. Col. (3) the 
observed 12 CO transition, Col. (4) the Gaussian beam FWHM, 
Col. (5) the lcr noise of the image. Col. ( 6 ) the conversion factor 
from flux density to brightness temperature for the three 12 CO 
transitions, and Col. (7) the /G 1 and R 32 line ratios. 


3.6. Star formation rate surface density maps 


Several SFR calibrations using different tracers, based on a vari¬ 
ety of galaxy an d/or region samples and stellar IMFs, have been 
publi s hed (e.g., Wu_eta]J l2005c ICalzetti et al.1 120071: IZhu et all 
120081 iRieke et al 20091) . As stated above, the Ha: recombina¬ 
tion emission line provides a nearly instantaneous measure of 
the SFR independently of the previous SF history. Among the 
hydrogen recombination lines, Her is the most widely used as 
SFR tracer because of its higher intensity and lower sensitivity 
to dust attenuation than bluer nebular lines (e.g., Fy a, Ly/i, H/j). 

To derive SFR maps we used the conversion between 
SFR and dust extinction-corrected H<r flux density derived by 
ICalzetti et al.l(l2007l) : 


2sfr = 5.3 x 10 - 42 S(Ha:) corr cos/ 


(5) 


where Xsfr is in units of M e yr 1 kpc “ 2 , and .S'(H<t) coit is the 
dust extinction-corrected Her in erg s -1 kpc -2 . This calibration 
has been derived assuming the stellar IMF of STARBURST99, 
which consists of two power laws with slope -1.3 in the range 
0.1 - 0.5 M G and slope -2.3 in the range 0.5 - 120 M 0 . (For de¬ 
tails on the adopted models of stellar populations see Appendix 
A 2 in ICalzetti et alJl2007l) 


We have also checked that the Her emission is uncontam¬ 
inated by emission from the AGN by comparing different (al¬ 
beit lower resolution) estimators of SFR. There is virtually no 
contamination from AGN-excited Ha' emission outside the cen¬ 
tral pixel for NGC 3627 and NGC 4826. Although NGC 4826 
has no X-ray source and NGC 3627 is only weakly detected 
(iHernandez-Garcla et alJ 120131) . NGC 4569 and NGC 4 579 
both have weak nuclear X-ray sources (iDudik et alJl2005h . In 
NGC 4579 potential AGN contamination may be more of a prob¬ 
lem because of the broad H a emission IIHo et al.lll997h . Never¬ 
theless, the SFR inferred from broad Ha' is ~0.02 M Q yr -1 , while 
the total SFR for this galaxy is ~2.2M e yr -1 (see Table0>, more 
than 100 times greater. We therefore conclude that SF processes 
are dominating the AGN in our sample galaxies. For safety, in 
NGC 4569 and NGC 4579 we masked the central 2x2 pixels 
since there was the possibility that the strong H» emission there 
is due to the AGN. 

3.7. Fitting method 

With the images described above, we constructed maps of Zh 2 
and Esfr to perform a pixel-by-pixel analysis of the spatially 
resolved K-S relation. We used data above 3cr significance both 
in Eh, and Z.sfr maps. 

We fit the data in logarithmic space: 

log(SsFR) = Ant + Nfit x log(E H ,) ( 6 ) 


where Afi t is the intercept and Nm the index of the K-S relation. 
We use the ordinary least squares (OLS) lin ear bisector method 


(jlsobe et al. 1990; Feigelson & Babu 

1992}), adopted in several 

SF relation studies (e.g., Bigiel et al. 

2008; Schruba et a 1.11201 It 

Momose et al. 2013|>. 


For statistical completeness, we also give results derived 
from a robust regression fitting method (l.i 1985c IFox II' l997), 
an alternative to least squares regression, which provides in¬ 
creased uncertainties of slopes and intercepts with respect to 
those given by the OFS bisector method. These fits are im- 
ple mented in the public-dom ain statistical software package 
R (llhaka & Gentlemanl 1 19961) . However, in the following, we 
present and discuss results emerging from the OFS bisector 
method for reasons described later in Sect. 16.31 

4. The molecular star formation relation for 
individual galaxies 

In this section we present the results obtained by the study of 
the pixel-by-pixel molecular K-S relation for each galaxy of our 
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Fig. 1: I H , contours in M G pc -2 derived from 12 CO(1-0), 12 CO(2-l), and 12 CO(3-2) overlaid on Zsfr images in M G yr 1 kpc -2 at 
the intrinsic spatial resolution of the 12 CO data and with pixel sizes equal to the 12 CO resolution. The beams are plotted (in yellow) 
in the lower left corners of maps and their values are listed in Table [2] (Acr, Ac))-offsets are with respect to the dynamical center of 
each galaxy. Eh, contours are drawn starting from noise levels >3cr with contour spacings that are whole multiples of >3cr. 
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sample, using the two or three lowest 12 CO transitions for deriv¬ 
ing E H , and Her emission for Esfr, with spatial resolution rang¬ 
ing from -20 pc to 200 pc. We show two sets of figures for each 
galaxy. The first one (Fig. Q} displays overlays of Eh 2 contours 
(above 3cr) in M Q pc -2 derived from 12 CO(1-0), 12 CO(2-l), and 
12 CO(3-2) emission line on Esfr maps in M e yr 1 kpc 2 derived 
from Her emission line at the resolution of original 12 CO data. 
The images involved in these overlays have been obtained fol¬ 
lowing prescriptions described in Sects. [3~4l 13.51 and 13.61 Each 
panel of Fig.|T]reports the galaxy name, the reference to the 12 CO 
emission line used to derive Eh, , and the spatial resolution under 
analysis, i.e. the one offered by the intrinsic 12 CO line map and 
at which the Her image has been degraded. 

The second set of images (Figs. 00 consists of two or three 
plots for galaxy, according to the available 12 CO line maps, dis¬ 
playing the K-S relation at the intrinsic resolution of the 12 CO 
line data. These plots contain the galaxy name, the reference 
to the adopted 12 CO line for deriving Eh,, the spatial resolu¬ 
tion, the fit line derived from the OLS bisector method de¬ 
scribed in Sect. 13.71 the values of the power index Wit and of 
the intercept Ag, and their uncertainties of the fit line accord¬ 
ing to equation 0, and the derived value for Td ep i defined as 
Tdepi - <Eh 2 )/<E sfr ) (e.g., [Leroy et al.|l2013l) . 

The available FoVs vary as a function of observed line and 
instrument. IRAM 12 CO(1-0) maps have a FoV of the primary 
beam of -42", IRAM 12 CO(2-l) of -21", SMA 12 CO(3-2) of 
-36", while HST Her images have usable FoVs ranging from 
-24" to -40". 

The results of this analysis, including the findings derived 
convolving all maps to the common spatial resolution of 200 pc, 
are collected in Table [4] In this table. Col. (1) indicates the 
galaxy name. Col. (2) the 12 CO transition used to derive Eh,, 
Col. (3) the spatial resolution in pc. Col. (4) the available FoV 
(diameter) in arcsec on the plane of the sky. Col. (5) the radius 
under investigation in kpc on the plane of the galaxy. Columns 
(6) and (7) give the power index ((Vg t ( OLSbis.)) and the intercept 
G4fit(OLS bis.)) and their uncertainties of the OLS bisector fitting 
line (for simplicity in the text these two parameters are indi¬ 
cated with Nn t and Ag t ). Column (8) indicates the Pearson cor¬ 
relation coefficient (r corr (OLSbis.)) of the OLS bisector fitting line 
(for simplicity in the text indicated with r co „) and the number of 
points under analysis (n. pts). Cols. (9) and (10) the power index 
(Vgt(RR)) and the intercept (AfjgRR)) and their uncertainties of the 
robust regression fitting line;. Col. (11) the mean Eh 2 (<Eh 2 >) 
in M 0 pc -2 within the FoV and taking only data points above 3<x 
significance into account. Col. (12) the mean Esfr (<Esfr>) in 
M 0 yr“ 1 kpc~ 2 under the same conditions. Col. (13) the molecu¬ 
lar Tdepi in Gyr, and Col. (14) the final pixel size of l2 CO and Ha 
images after procedures described in Sect. 13.41 

4.1. NGC 3627 

NGC 3627 (Messier 66) is an interacting (e.g., ICasasola et al.l 
120041) and barred galaxy classified as SAB(s)b at a distance 
of 10 Mpc, with signatures of a LINER/Seyfert 2 type nu¬ 
clear activity dHo et al.il9971) . With NGC 3623 and NGC 3628, 
it forms the well-known Leo Triplet (M 66 Group, VV 308). 
Optical broad-band images of NGC 3627 reveal a pronounced 
and asymmetric spiral pattern with heavy dus_t_lanes, indicat¬ 
ing strong density wave action dPtak et al.l 120061) . While the 
western arm is accompanied by weak traces of SF visible in 
Ha, the e astern arm contain s a star-forming segm ent in its in¬ 
ner part dSmith et al.l Il994t IChemin et al.l l2003h . NGC 3627 
also possesses X-ray properties of a galaxy with a recent star- 


burst (Dahlemetal. 1996). Both the radio continuum (I 

1.8 cm 
nd the 

and 20 cm, Urbanik et al. 1985; 

Paladino et al. 2008) a 

CO emissions (e.g., Regan et al 

2001; Paladino et al. 

2008| 

Haan et al. 2009; Casasola et al. 

2011; Watanabe et al. 

2011b 

show a nuclear peak, extending along the leading edges of the 


bar that forms two broad maxima at the bar ends, and then the 
spiral arms trail off from the bar ends. In contrast, the Hi emis¬ 
sion exhibits a s piral morphology without signatures o f a bar in 
the atomic gas (lHaan et ail 120081 : 1 Walter et al.ll2008l) . The de¬ 
rived gravity torque budget shows that NGC 3627 is a poten- 
tialjimoking gun of inner gas inflow at a resolution of -60 pc 
(ICasasola et al J201lh . In addition to 12 CO lines, other molecular 
transitions have been detected in NGC 3627, including HCN(1- 
0), HCN(2-1), HCN(3-2), HCO + (1-0), and HCO + (3-2), sug¬ 
gesting the presence of high density gas (e.g., iGao & Solomonl 
l2fi04tlKrins etakll2008lf . 

The panels of the first line of Fig. Q] show the superposition 
of Eh, image contours derived from 12 CO(1-0) and 12 CO(2- 
1) emission lines overlaid on Esfr images estimated from the 
Ha emission at the intrinsic resolutions of CO maps, i.e., —1'.'7 
(-82 pc) for the (1-0) transition and ~0'/7 (-36 pc) for the (2- 
1) one. The value of Eh, has been derived from 12 CO(1-0) and 
12 CO(2-l) emission lines by using Eqs. (0 and 0, respectively, 
and Esfr by using Eq. 0. For l2 CO(l-0) (left panel), the Eh 2 
and Esfr peaks are spatially coincident in the limit of the res¬ 
olution of —1"7 and the two distributions are quite consistent 
within a radius of -3—4" (150-200pc) from the nucleus (on the 
plane of the sky). At larger distances. Eh, and Esfr are not cor¬ 
related, mainly because Esfr is not distributed along a bar as is 
Eh, • Similar (anti-)correlations also characterize the comparison 
between Eh, derived from 12 CO(2-l) emission line and Esfr at 
the resolution of -077 (right panel). 

The lefthand panel of Fig. [2] shows the molecular K-S rela¬ 
tion derived for NGC 3627 using the 12 CO(1-0) emission line 
to estimate Eh, at the resolution of 82 pc. At this resolution the 
power index Vg t is equal to 1.18 + 0.11 within a radius of 1.3 kpc 
(on the plane of the galaxy). Following the same procedure but 
convolving the original l2 CO(l-0) and Ha maps at the lower res¬ 
olution of 200 pc, the K-S relation has A/,,, equal to 1.11 ± 0.41. 
The righthand panel of Fig. [2] shows the results obtained using 
the 12 CO(2-l) emission line to derive En,(Eq.0 at the resolu¬ 
tion of 36 pc. At this resolution /V(; t = 1.16 + 0.05 within a radius 
of 1.1 kpc (on the plane of the galaxy), while /V (l , increases to 
1.59 ± 0.63 at the resolution of 200 pc. Although all the derived 
/Vf,| values are consistent with literature results (see later discus¬ 
sion in Sect. ED- the K-S relations studied at the resolution of 
200 pc - for both 12 CO(1-0) and 12 CO(2-l) emission lines - 
with only four data points involved in the analysis do not allow 
us to infer statistical conclusions. In NGC 3627, we can only say 
that Vet « 1.2 both with Eh, derived from 12 CO(1-0) at 82 pc 
and from 12 CO(2-l) at 36 pc. 

Neglecting the 200 pc-resolution cases, the Pearson correla¬ 
tion coefficient is -0.6 - 0.7 and Td ep i is -1.2 - 1.3 Gyr for both 
the lowest 12 CO transitions. 

4.2. NGC 4569 

NGC 4569 (Messier 90) is a bright SAB(rs)ab galaxy at a dis¬ 
tance of 17 M pc in the Virgo cluster. A l arge scale bar is seen 
in NIR images dLaurikainen & Sakll2002l) and is almost aligned 
with the maj or axis of the galaxy (PA = 15 deg according to 
Jo gee et a l.l2005l). The galaxy ha rbors a nucleus of the transition 
type (type T2 in iHo et all 1 1 997l) . which exhibits a pronounced 
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Fig. 2: Left panel: The K-S relation plot for NGC 3627 at the resolution of 82 pc. E H , was derived from the l2 CO(l-0) emission line 
map and Esfr from the Her image based on Eqs. Q and (0, respectively. Red triangles indicate data points above 3<x significance 
both in E H , and Esfr, within a radius of 1.3 kpc (on the plane of the galaxy). The solid black line indicates the OLS bisector fitting 
line (see Sect. I3.7l for fitting method). The index A% t and the intercept ,4n, of the OLS bisector fitting line, the Pearson correlation 
coefficient r core , and r,| L .pi values are reported in figure. Right panel: Same as left panel with In, derived from the 12 CO(2-l) emission 
line map based on Eq. 0 at the resolution of 36 pc. Red triangles indicate data points above 3cr significance both in Eh, and Esfr, 
within a radius of 1.1 kpc (on the plane of the galaxy). 




Fig. 3: Left panel: Same as left panel of Fig.[2]for NGC 4569 at the resolution of 154 pc. Red triangles indicate data points above 3cr 
significance both in Eh, and Esfr, within a radius of 2.9 kpc (on the plane of the galaxy). Right panel: Same as left panel with E H , 
derived from the 12 CO(2-l) emission line map at the resolution of 74 pc. Red triangles indicate data points above 3<x significance 
both in Eh, and Esfr, within a radius of 2.1 kpc (on the plane of the galaxy). 


nuclear starburst activity. Hi emission line observations have re¬ 
vealed that in NGC 4569 the atomic gas is distributed along a 
central bar with radius of -60" (~5 knc. lHaan et al.l2008h . Inter- 
fero metric CO observations of NGC 4569 presented previously 

(e.gJHelfer et al. 2003lljogee et all2005UNakanishi et all2005t 

Boo ne et al.ll2007 ) have shown that the major part of the molec¬ 
ular gas detected in the inner 20" is concentrated within a radius 
of 800 pc, distributed along the large scale stellar bar seen in NIR 
observations, and with a peak close to the center and another 
one at —500 pc from it. A hole in the CO distribution coincides 
with the nu cleus where most of the Her emission a nd blue light 
are emitted dPogge et al.ll2000l) . lBoone et al.l (120071) also demon¬ 


strate that the gravitational torques are able to efficiently funnel 
the gas down to -300 pc. 

The panels of the second line of Fig. Q] show the compari¬ 
son between Eh. and Esfr distributions in NGC 4569 at resolu¬ 
tions of —1'/9 (-154 pc from 12 CO(1-0)), -079 (-74 pc from 
12 CO(2-l)), and ~2'.'3 (-190 pc from 12 CO(3-2)). The case in¬ 
volving the 12 CO(3-2) emission line refers to Eq. © for deriving 
Eh 2 . >From these overlays it can be seen that Eh, and Esfr are 
differently distributed at the available resolutions. While Esfr 
is centrally concentrated and peaked. Eh, (derived both from 
12 CO(1-0) and l2 CO(2-l) emission line) is distributed along a 
large scale bar of -17" x 6" (-1.4 kpc x 0.5 kpc) in size (in 
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the plane of the sky) whose peak is far away from the peak of 
Isfr. The same distribution is also visible in the SMA 12 CO(3- 
2) map. There are, however, differences, perhaps due to the lower 
S/N ratio of the 12 CO(3-2) SMA data relative to the IRAM ob¬ 
servations (iBoone et al.ll201 lb . 

Similar to NGC 3627, NGC 4569 exhibits a centrally con¬ 
centrated morphology in the E$fr images and a bar-like dis¬ 
tribution in the E gas maps. This could mean that the action of 
the bar has transported the gas to the nuclear regions to fuel a 
mini-starburst episode as found in many nearby spiral galaxies 
dSakamoto et al.fi 999)) . 

Figure [3] and the lefthand panel of Fig.[6]show the molecular 
K-S relation derived for NGC 4569 with Eh, estimated from the 
three lowest 12 CO emission lines at the intrinsic resolution of 
the 12 CO maps. These figures, together with findings collected 
in Table [4] show that NGC 4569 has an index AW of the K-S 
relation that is sublinear (-0.6 - 0.7) for the three available 12 CO 
transitions studied at spatial scales from 74 to 200 pc (producing 
six subcases), roughly constant as a function of resolution for a 
given l2 CO transition, and only slightly varying as a function of 
12 CO transition at a resolution of 200 pc. 

For NGC 4569, the Pearson correlation coefficient is ap¬ 
proximately invariant with respect to the resolution, from 74 to 
200 pc, for a given l2 CO transition. In contrast to this, r corr varies 
as a function of l2 CO line at a given resolution (i.e., 200 pc). The 
best r corr is obtained with the 12 CO(2-l) line (-0.5-0.6), while 
the worst one with the 12 CO(3-2) transition (-0.3) ( 12 CO(1-0) 
gives r mir -0.4). 

As in NGC 3627, NGC 4569 has a short molecular Td ep i of 
-1 Gyr, suggesting that the gas is efficiently converted in stars. 
Within a radius of 0.5 kpc, Td ep i is even smaller assuming values 
of -0.7-0.9 Gyr. 


4.3. NGC 4579 


NGC 4579 (Messier 58) is a SAB(rs)b galaxy classified as an 
intermediate type 1 object (LINER/Seyfert 1.9) by iHcTetal] 
d 19971) at a distance of 20 Mpc in the Virgo cluster. It 
also has an unresolved nuclear hard X-ray (variable) source 
with a prominent broad Fe Ka line dTerashimaetalJj2000; 
iHo et all 1200 lb lEracleous et alJ l2002t iDewangan et al] 12004). 
A non-thermal radi o continuum source is detected at the po¬ 
sition of the AGN ( Hummel etal. 19871 : IHo & Ulvestadll200lb 
lUlvestad & Hol2001blKrips et al.ll2007l) . The NIR K-band image 
of NGC 457 9 has revealed a large-scal e stellar bar and a weak 
nuclear oval dGarcla-Burillo et al]j2009i) . The 21 cm Hi line ob¬ 
servations have shown that the atomic gas in NGC 4579 is cur¬ 
rently piling up in a pseudo-ring (radius of -40", -4 kpc) formed 
by two winding spiral arms that are morphologically decoupled 
from the bar structure ( lHaan et al]|2008l) . Molecular gas in the 
inner r < 2 kpc disk is distributed in two spiral arms , an outer arc 
and a central lopsided disk-like structure (lGarc fa-Burill o~et~akl 
l2009t) . The derived gravity torque budget in NGC 4579 have 
shown that inward gas flow is occurring on different spatial 
scales in the disk, with clear smoking gun evidence of inward 
gas transport down to r - 50 pc. 

The panels of the third line of Fig.Q]show that in NGC 4579, 
Eh, and Esfr have completely different distributions distribu¬ 
tions at resolutions of both ~1'.'6 (-154 pc from 12 CO(1-0)) 
and -078 (-76 pc from 12 CO(2-l)). The distribution of Esfr 
is centrally concentrated within the inner -4" (on the plane of 
the sky) and, at larger distances from the center, lies along two 
spiral arms extending up to -12" from the nucleus. The mor¬ 
phology of E h , traced by 12 CO(1-0) and 12 CO(2-l) is instead 


mainly defin ed by two highly contraste d spiral lanes without a 
central peak dGarcla-Burillo et ahll2009l) . 

Figure [4] shows the results obtained from the analysis of the 
K-S relation for NGC 4579. The value of AW ranges from -0.5 
to -1.1 on spatial scales of 76-200 pc and 

considering 12 CO(1-0) and 12 CO(2-l) lines for E H , deriva¬ 
tion. For a given 12 CO line, Vat decreases with finer resolu¬ 
tion. While for l2 CO(l-0) Wit gradually decreases from -0.9 
at 200 pc-resolution to -0.8 at 154 pc-resolution, for 12 CO(2-l) 
the decreasing of Wit is stronger, from 1.1 to 0.5, possibly be¬ 
cause the change from 200 to 76 pc in resolution is more drastic. 
At the common resolution of 200 pc, AW increases with higher 
/-CO transition. 

In NGC 4579, the quality of the correlation tends to worsen 
with resolution. While for l2 CO(l-0) r corr is approximately con¬ 
stant (-0.6) through 154-200 pc resolution, for 12 CO(2-l) it 
drops down from 0.6 to 0.5 in the resolution range from 200 to 
76 pc. The correlation coefficient r corr does not seem to depend 
on the CO transition, since it is -0.6 for both 12 CO lines at the 
common resolution of 200 pc. 

Unlike NGC 3627 and NGC 4569, NGC 4579 has a more 
“standard” molecular Td ep i of -2 Gyr (see later discussion in 
Sect. 15.21) . without a radial trend. The similar trends obtained for 
AW and r con as a function of the resolution and for both the two 
lowest l2 CO transitions suggest that in NGC 4579 the K-S rela- 
tionis almost invariant with respect to 12 CO(1-0) and 12 CO(2-l) 
emission lines, but it changes slightly with spatial resolution. 

4.4. NGC 4826 

NGC 4826, also known as the “Black Eye” or “Evil Eye” galaxy 
due to its optical appearance, is the closest target (-5 Mpc) of the 
NUGA core sample, and its nucleus is classified as a LINER type 
(IHo et al 111 9971) . It hosts two nested counter-rotatin g atomic and 
molecular gas disks of comparable mass (-10 8 M„. iBraun et alJ 
1 1992b ICasoli & Gerinlfl993t IBraun et aD 1 19941) . The inner disk 
has a radius of ~50"(~1.3 kpc), while the outer one from -80" 
to -9'8 (-2.1 - 15.3 kpc). By studying the_stellar kinematics 
along the principal axes of NGC 4826. iRix et al.l (1 19951) found 
that the stars rotate at all radii in the same direction as the inner 
disk, providing strong evidence that stars and gas are coplanar. 
NUGA observations have shown a high concentration of molec¬ 
ular gas (in 12 CO(1-0), 12 CO(2-l), and 12 CO(3-2)), within a 
radius of 80 pc, forming a circumnuclear molecular disk. A de¬ 
tailed analysis of the kinematics, however, does not reveal any 
evidence of fueling of the nucleus. 

The panels of the fourth line of Fig. Q] show the comparison 
between Eh, and Esfr distributions in NGC 4826 at a resolution 
of -271 (-56 pc from 12 CO(1-0)), -077 (-17 pc from 12 CO(2- 
1)), and -273 (-59 pc from 12 CO(3-2)). Unlike NGC 4569 and 
NGC 4579, the morphologies of Esfr and Eh, in NGC 4826 are 
spatially coincident and characterized by a structured disk. 

Figure[5]and the righthand panel of Fig.[6]show the molecu¬ 
lar K-S relation derived for NGC 4826 with Eh, estimated from 
the three lowest 12 CO emission lines at the intrinsic resolution of 
the gas maps. As in NGC 3627, there are some cases where only 
four or five data points are involved in the analysis ( 12 CO(2-l) 
and 12 CO(3-2) at 200 pc). Neglecting these cases, the K-S re¬ 
lation derived from 12 CO(1-0) shows that AW decreases from 
-1.3 at 200 pc-resolution to -0.9 at 56 pc-resolution. The case 
with 12 CO(1-0) at 200 pc has a superlinear AW, while those an¬ 
alyzed at the intrinsic resolution of the all three 12 CO maps (17- 
59 pc) are sublinear. The comparison between the K-S relation 
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Fig. 4: Left panel: Same as left panel of Fig.[2]for NGC 4579 at the resolution of 154 pc. Red triangles indicate data points above 3cr 
significance both in Eh, and Esfr, within a radius of 2.7 kpc (on the plane of the galaxy). Right panel: Same as left panel with X H , 
derived from the 12 CO(2-l) emission line map at the resolution of 76 pc. Red triangles indicate data points above 3cr significance 
both in Eh, and Esfr, within a radius of 1.5 kpc (on the plane of the galaxy). 




Fig. 5: Left panel: Same as left panel of Fig.[2]for NGC 4826 at the resolution of 56 pc. Red triangles indicate data points above 3<x 
significance both in Eh, and Esfr, within a radius of 1.1 kpc (on the plane of the galaxy). Right panel: Same as left panel with Eh 2 
derived from the 12 CO(2-l) emission line map at the resolution of 17 pc. Red triangles indicate data points above 3<x significance 
both in Eh, and Esfr, within a radius of 0.6 kpc (on the plane of the galaxy). 


with 12 CO(1-0) at 56 pc and that with 12 CO(3-2) at 59 pc gives 
an almost identical tVfi t (—0.9). 

The good agreement between Esfr and Eh, distributions 
(Fig. ID consequently corresponds to a high Pearson correlation 
coefficient, equal to ~0.7-0.8 for all the cases with more than 
five data points at resolutions of 56-200 pc. These values of r con 
are the best in our sample of galaxies, and those of r corr drop 
down to ~0.5 in the case using the 12 CO(2-l) line at 17 pc, the 
highest resolution available not only in NGC 4826 but also in the 
present study. This drop of r corr could be due to the high resolu¬ 
tion; however, the comparison with the K-S relationship derived 
using the 12 CO(2-l) line at 200 pc is not definitive proof of this 
because of the poor number statistics. For other sample galaxies 
(NGC 4569 and NGC 4579), there is no a strong difference in 
r c(m between l2 CO(l-0) and 12 CO(2-l). In addition to A^ t , r co „ 


also assumes similar values (0.7-0.8) in the cases of l2 CO(l-0) 
at 56 pc and 12 CO(3-2) at 59 pc. The results obtained for ;V (lt 
and r con therefore suggest that NGC 4826 has a K-S relation¬ 
ship independent of the 12 CO transition at resolutions as high as 
~60 pc. 


Like NGC 4579, NGC 4826 has a Td ep i of ~2 Gyr without 
radial trend both by using 12 CO(1-0) and 12 CO(2-l), when tak¬ 
ing only the statistically significant cases into account. The use 
of the 12 CO(3-2) line instead gives a shorter Tdepi, of ~ 1 Gyr, 
and radial trend since it goes from ~ 1 Gyr within r < 1.0 kpc to 
~0.7 Gyr within r < 0.2 kpc. 
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Table 4: Derived parameters from the pixel-by-pixel molecular star-formation relation for each galaxy. 


Galaxy 

(1) 

lz CO transition 

(2) 

Spatial resolution 
[pc] 

(3) 

FoV sky 

["] 

(4) 

fgal 

[kpc] 

(5) 

Vfit(OLSbis-) 

(6) 

Afit(OLSbis.) 

(7) 

rcorr(OLS bis.) (n. pts) 

(8) 

-VfitfRR) 

(9) 

Afit(RR) 

(10) 

<E H 2> 
[Mo pc~ 2 ] 
(11) 

<^SFR> 

[Moyr-'kptT 2 ] 

(12) 

Tdepl 

[Gyr] 

(13) 

scale 

r/pix] 

(14) 

NGC 3627 

(1-0) 

82 

24 

1.3 

1.18 ±0.11 

-3.66 ±0.19 

0.69 (12) 

0.83 ±0.27 

-2.72 ± 0.73 

50.69 

0.04 

1.17 

2.0 


(1-0) 

200 

24 

1.3 

1.11 ±0.41 

-3.47 ± 0.74 

0.70 (4) 

0.75 ±0.18 

-2.47 ± 0.47 

53.60 

0.04 

1.49 

4.0 


(2-1) 

36 

20 

1.1 

1.16 ±0.05 

-3.70 ± 0.07 

0.62 (88) 

0.72 ±0.10 

-2.63 ±0.27 

82.32 

0.06 

1.27 

0.8 


(2-1) 

200 

20 

1.1 

1.59 ±0.63 

-4.76 ± 0.92 

0.65 (4) 

1.06 ± 1.05 

-3.32 ±2.79 

75.19 

0.05 

1.45 

4.0 

NGC 4569 

(1-0) 

154 

30 

3.1 

0.60 ± 0.03 

-1.96 ±0.03 

0.42 (37) 

0.18 ±0.08 

-1.02 ±0.19 

51.40 

0.16 

1.03 

2.0 


(1-0) 

200 

30 

3.1 

0.61 ± 0.05 

-2.02 ± 0.06 

0.43 (30) 

0.24 ±0.10 

-1.21 ±0.22 

50.43 

0.21 

0.99 

2.5 


(2-1) 

74 

20 

2.1 

0.59 ± 0.01 

-1.99 ±0.01 

0.59(145) 

0.32 ± 0.04 

-1.37 ±0.09 

116.29 

0.36 

1.17 

1.0 


(2-1) 

200 

20 

2.1 

0.64 ± 0.05 

-2.11 ±0.11 

0.54 (20) 

0.35 ±0.13 

-1.45 ±0.30 

92.68 

0.31 

1.12 

2.5 


(3-2) 

190 

30 

2.9 

0.67 ± 0.04 

-2.23 ±0.11 

0.27 (39) 

0.07 ± 0.06 

-0.96 ±0.14 

46.50 

0.27 

1.05 

2.5 


(3-2) 

200 

30 

2.9 

0.69 ± 0.08 

-2.28 ± 0.08 

0.30 (36) 

0.08 ± 0.07 

-0.99 ±0.16 

42.63 

0.25 

1.03 

2.5 

NGC 4579 

(1-0) 

154 

36 

2.7 

0.80 ± 0.04 

-2.94 ± 0.02 

0.60 (65) 

0.47 ± 0.08 

-2.27 ±0.16 

17.36 

0.01 

2.28 

1.8 


(1-0) 

200 

36 

2.7 

0.91 ±0.01 

-3.13 ±0.02 

0.58 (61) 

0.57 ± 0.09 

-2.48 ±0.18 

17.05 

0.01 

2.07 

2.0 


(2-1) 

76 

20 

1.5 

0.54 ± 0.01 

-2.29 ± 0.02 

0.49(111) 

0.24 ± 0.04 

-1.69 ±0.08 

40.64 

0.02 

2.22 

1.0 


(2-1) 

200 

20 

1.5 

1.14 ±0.03 

-3.49 ± 0.03 

0.62 (40) 

0.71 ±0.16 

-2.66 ±0.31 

37.76 

0.02 

1.60 

2.0 

NGC 4826 

(1-0) 

56 

40 

1.1 

0.87 ± 0.03 

-2.94 ± 0.03 

0.72 (70) 

0.64 ± 0.07 

-2.43 ±0.17 

41.38 

0.02 

1.77 

2.0 


(1-0) 

200 

40**9 

1.1 

1.27 ±0.08 

-3.89 ± 0.05 

0.77 (12) 

0.91 ±0.25 

-3.09 ±0.53 

105.07 

0.05 

1.98 

7.6 


(2-1) 

17 

20 

0.6 

0.67 ± 0.01 

-2.45 ± 0.01 

0.52 (158) 

0.33 ±0.04 

-1.71 ±0.10 

112.31 

0.05 

2.07 

1.0 


(2-1) 

200 

20 (b) 

0.6 

1.64 ±1.31 

-4.17 ±0.86 

0.49 (4) 

0.90 ± 1.13 

-2.71 ±2.26 

99.85 

0.13 

0.79 

7.6 


(3-2) 

59 

36 

1.0 

0.89 ± 0.04 

-2.73 ± 0.07 

0.76 (27) 

0.64 ±0.11 

-2.19 ±0.24 

16.28 

0.02 

1.07 

2.0 


(3-2) 

200 

36 (c) 

1.0 

1.22 ±0.11 

-3.19 ±0.17 

0.91 (5) 

1.11 ±0.30 

-3.00 ±0.53 

20.93 

0.03 

0.61 

7.6 


Notes. ( " ) The effective FoV is reduced to 30'.'4 since we rebinned 40" with a pixel scale of 776. tb) The effective FoV is reduced to 1572 since we rebinned 20" with a pixel scale of 776. <rl The 
effective FoV is reduced to 3074 since we rebinned 36" with a pixel scale of 776. 
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Fig. 6: Left panel. Same as Fig. [3] for NGC 4569, with Ih 2 derived from the 12 CO(3-2) emission line map based on equation (4} at 
the resolution of 190 pc. Red triangles indicate data points above 3cr significance both in I H , and Isfr, within a radius of 2.9 kpc 
(on the plane of the galaxy). Right panel: Same as left panel for NGC 4826 at the resolution of 59 pc. Red triangles indicate data 
points above 3cr significance both in Ih 2 and Is hr, within a radius of 1.0 kpc (on the plane of the galaxy). 


Table 5: Derived parameters from the pixel-by-pixel molecular 
star-formation relation across the sample at the common spatial 
resolution of 200 pc. 


12 CO transition 

( 1 ) 


Nat 

( 2 ) 


A m 
(3) 


r (n.pts) 

(4) 


( 1 - 0 ) 


1.20 + 0.03 -3.60 + 0.01 0.66(107) 


( 2 - 1 ) 


1.14 + 0.01 -3.41 +0.02 0.77 (68) 


(3-2) 


0.67 ± 0.02 -2.22 ± 0.08 0.41 (41) 


(l-0)+(2-l)+(3-2) 1.14 + 0.01 -3.41+0.01 0.64(216) 


5. The molecular star-formation relation across the 
sample 

The analysis performed for each galaxy showed a wide range 
of behaviors in terms of the K-S relations, and there is no “uni¬ 
versal” molecular SF relation although the investigated galaxies 
belong to the same subclass of objects (i.e., nearby active galax¬ 
ies), and all the derived quantities have been treated with the 
same methodology and at comparable spatial resolutions. The 
main result is therefore that each galaxy has its own SF rela¬ 
tion (with its own index Nr t , correlation coefficient, and Td ep i) on 
spatial scales of -20-200 pc. Nevertheless, we identified some 
common behaviors in terms of K-S relation, as discussed below. 

5.1. The K-S relation index: resolution vs. 12 CO transitions 

By using the three lowest I2 CO at resolutions of -20-200 pc, 
we found K-S relation indexes Ajj, ranging from —0.5 to - 1.6, 
which are all values that are consistent with literature results. A 
superlinear slope of the K-S relation is consistent both with the 
earl y global di s k stu dies of the SF relation bv lKennicuttl (fl998a) 
and IKennicutt! (1 1 998hl) based on the combination of atomic and 
molecular gas data for the I gas computation and with more re¬ 
cent works based only on molecular component at sub-kpc scales 
(e.g., iKennicutt et all 120071 IVerlev et all 120101: ILiu et alJ 1201 lb 


iRahman et alJ 1201 it iMomose et alJ l2013h . iBigiel et aD (2008) 
instead derived N « 1 from the correlation between Isfr and Ih 2 
estimated from 12 CO(2-l) data for seven nearby spiral galaxies. 
As already mentioned in Sect.lTl lBigiel et akld2008l) suggest that 
_a linear correlation is evident in regions of high gas surface den¬ 
sities where the gas is typically molecular (>10 M Q pc~ 2 ). Other 
studies using 12 CO(2-l) also show a linear correlation (e.g., 
iLerov et al .1120081 ISchruba et al.ll20l71) . including one that com- 
""bined si ngle-dish 12 CO (2-l ) data with interferometric 12 CO(l- 
_ 0) data (IRahman et a 1 .11 20111) . These 12 CO(2-l) studies analyzed 
_ a substantial number of nearby galaxies, though it should be rec¬ 
ognized that some studies based on 12 CO(1-0) data showed a 
superlinea r (power-law) correlation, rather than a linear correla¬ 
tion ( e.g., IWong & Blitzl 12002k IKennicutt et alJl2007t ILiu et alJ 
1201 lh . This suggests that the choice of the CO transition for 
deriving In, could affect Ag t valu es (see Bjgigietal. 2008). A 
linear slope has also been found bv iVutisalchavakul et al.f (i2014t) 
for the SF relation in an 11 deg 2 region of the Galactic plane with 
dust continua at 1.1 mm and 22 pm emission used as tracers of 
molecular gas and SFR, respectively, over a range of resolution 
from 33" to 20'(-0.1-45 pc). 

The most common explanation for a linear K-S relationship 
is that the observed CO luminosity is directly proportional to the 
number of star-forming clouds or GMCs, with all clouds having 
similar properties, such as the volume density, the efficiency of 
the cloud, and the SFR. In observations at resolution k. 100 pc, 
the individual clouds are not resolved but rather their CO flux is 
dispersed throughout the beam. In this case, regions with more 
clouds emit more CO in proportion to the number of clouds. 
Other recent works favor a subline ar (N - 0.6-0.8) K-S relation¬ 
ship at resolutions > 170 pc (e .g., iBlanc et akll20()9 j : iFord et all 
120 1 3t Ishettv et al,ll2013ll2014l) , A sublinear K-S relationship, in 
contrast to a linear one, suggests that the clouds do not have 
the same properties, and SFRs and/or volume densities vary 
( iShettv et al.l2014l) . This means that there is no one-to-one corre¬ 
spondence between the CO luminosity and the number of clouds. 
In addition, the conversion factor .Xco also varies with the loca¬ 
tion in a galaxy (see later Sect. ED- Finally, a more straight¬ 
forward explanation for the sublinear K-S relationship would be 
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Fig. 7: K-S relation plot for all sample galaxies. The top left panel plots all galaxies and only 12 CO(1-0) emission line data; the 
top right panel all galaxies and only 12 CO(2-l) emission line data; the bottom left panel all galaxies and only 12 CO(3-2) emission 
line data; and the bottom right panel all galaxies and all the available 12 CO emission line data, at the common spatial resolution of 
200 pc. Different symbols and/or colors indicate galaxies whose Eh, has been derived from a given l2 CO transition, as shown in 
the labels. The black lines refer to the OLS fits derived for different cases: the short-dotted line is for all galaxies and l2 CO(l-0) 
emission line data, the point-dotted line for all galaxies and 12 CO(2-l), the long-dotted line for all galaxies and 12 CO(3-2), and the 
black solid line is for all galaxies and all the available 12 CO lines. In the bottom right panel the 12 CO(2-l) OLS fit line is hidden by 
total OLS fit line. 


that CO permeates the hierarchical interstellar medium, includ¬ 
ing the filamen ts and lower density regi ons within which GMCs 
are embedded (IShettv et alJl2013l.l2014l) . 

>From Table [4] it emerges that for a given l2 CO transition, 
the index A/u tends to gradually decrease with finer resolution 
down to a resolution of approximately 20 pc. On the other hand, 
on scales larger than ~100 pc, the slope of the SF relation tends 
to increase and become somewhat superlinear. This result is 
common to all three lowest 12 CO transitions and to all galaxies 
except for NGC 4569 whose A/a t is almost constant as a func¬ 
tion of resolution as shown in the previous section. Our study 
therefore suggests that ~80-100 pc is the scale at which the K-S 
relation undergoes a change, both with respect to kpc scales and 
slightly finer spatial scales (<250 pc, e.g jBigiel et alj '2008. and 
this work). This contrasts with previous results by Onodera et al.l 
( 2010 ) for the Local Group spiral galaxy M 33, in which the 


slope steepens, and the correlation virtually disappears on small 
spatial scales of ~80 pc (i.e., GMC scales). However, in M 33, 
both the SFR and gas surface densities at 80 pc resolution are 
significantly lower than in the NUGA sample. 

For distinguishing effects on the K-S relation index of dif¬ 
ferent 12 CO transitions to derive Eh, , data are needed at a com¬ 
mon resolution, 200 pc in the present work. The results from this 
comparison are displayed and collected in Fig. |7] and Table [3 
respectively. Each panel of Fig. |7] shows the K-S relation plot 
for all sample galaxies and taking the three 12 CO lines into ac¬ 
count both separately and all together without distinction based 
on the transition. Different symbols and/or colors indicate galax¬ 
ies whose Eh, has been derived from a given 12 CO transition, 
and the black lines refer to the OLS fits derived for different 
cases. At the resolution of 200 pc, taking all sample galaxies 
into account and 12 CO(1-0) (top left panel) and l2 CO(2-l) (top 
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right panel) separately, the resulting K-S relations have indices 
of 1.20 + 0.03 and 1.14 + 0.01, respectively. These findings in¬ 
dicate that 12 CO(1-0) and 12 CO(2-l) lines lead to quite similar 
results in terms of the slope of the K-S relation. Under the same 
conditions but by using 12 CO(3-2) data (bottom left panel), we 
found a sublinear slope of 0.67 + 0.02. However, this result is 
based only on two cases, NGC 4569 and NGC 4826, and for 
NGC 4569 all slopes for the three lowest 12 CO transitions at 
all investigated resolutions (74-200 pc) are sublinear. In con¬ 
trast, in NGC 4826 the slope derived from the 12 CO(3-2) line is 
roughly unity at 200 pc of resolution (see Table HJ). Taking all 
sample galaxies and all available 12 CO transitions (bottom left 
panel) into account, the K-S relation has a “standard” index of 
1.14 + 0.01. 

It is well known that each transition traces different physical 
gas properties. The kinematic tempe rature of molecular gas is 
typically ~10 K (IScoville et alJl 1987h . which is above the level 
energy temperature of 5.5 K for the / = 1 level of the 12 CO, 
but below the temperatures of 16.5 K and 33 K for the J = 2 
and J - 3 levels, respectively. This implies that a slight change 
in gas kinematic temperature is sufficient to affect the excita¬ 
tion for the 12 CO(2-l) and 12 CO(3-2) emission lines. The dif¬ 
ferent critical densities of the three 12 CO transitions (~10 3 cm -3 , 
~2 x 10 4 cnr 3 , and ~7 x 10 4 cnr 3 for 12 CO(1-0), 12 CO(2- 
1), and 12 CO(3-2), respectively) make their line ratios sensi¬ 
tive to local gas density. The value of /Gi has been observed 
to systematically vary with SFR and I gas in the Milky Way_(e.g., 
Sawada et al. 1200 li) and M 51 (Koda et ali 
201 3). Despite this, because of the use 
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of individually measured mean Kyi to convert to an equivalent 
CO(l-O) mass conversion, such variations do not have a signifi¬ 
cant impact on the resulting index N&t of the K-S relation. 

Another result emerging from Table |4]concerns the K-S rela¬ 
tion studied at 200 pc resolution with I H2 derived from 12 CO(l- 
0). NGC 3627, NGC 4579, and NGC 4826 exhibit K-S relations 
with indices /V (l , around unity, with In, of ~20-100M o pc~ 2 , and 
+depi of ~ 1.5-2 Gyr. NGC 4569 is different from the other sample 
galaxies also in this respect; it indeed has a I H , of -50 M 0 pc -2 , 
comparable to the other galaxies, but its K-S relation has a sub¬ 
linear slope and Td ep i of ~1 Gyr. 

5.2. The molecular gas depletion time 

As discussed in Sect. Q] r ( i ep i is an important parameter in the 
study of the K-S relation, as is the power index N& t , for under¬ 
standing the SF phenomenon in galaxies. We found molecular 
Tdepi ranging from ~1 to ~2 Gyr, which is consistent with more 
recent results as discussed below. Each galaxy has its own Td ep i 
that seems to be invariant with respect to the spatial scale probed 
in a given galaxy, at resolutions from ~20 to 200 pc. The molecu¬ 
lar Tdepi is also invariant with respect to 12 CO(1-0) and 12 CO(2- 
1) lines. The discussion of the Td ep i derived from 12 CO(3-2), for 
NGC 4569 and NGC 4826, deserves separate treatment and is 
discussed later. 

Both NGC 3627 and NGC 4569 have molecular Td ep i of ~1- 
1.3 Gyr, by using 12 CO(1-0) and 12 CO(2-l) and taking only the 
statistically significant cases into account. These values of Td ep i 
are compatible with the mean molecular [ I 2 CO( 1—0)] Td ep i of 
~1.2 Gyr (with aco = 3.5 M 0 (K km s~' pc~ 2 )~') found by 
ISaintonge et al.l ( 2011 ) for the COLD GASS sample of more 
than 200 galaxies at distances of -100-200 Mpc, but lower 
than the me a n mol ecula r f 12 CO(2-l ) l Td Rn i of ~2 Gyr found by 
Bigie l et al.l (120081) and iLerov et al.l (f2008l) in the HERACLES 


survey consisting of a sample of 48 nearby (D ~ 3-20 Mpc) 
galaxies. NGC 3627 and NGC 4569 therefore have a Td ep i that 
is more consistent with the COLD GASS single-dish survey that 
gives a global picture of gas and SF in the Local Universe but 
lacks the power to trace the exact distribution of these com¬ 
ponents, rather than with the HERACLES survey, which like 
the present study, maps them. However. ISaintonge et al.l ( 2011 ) 
demonstrate that after using the same conversion factor from CO 
luminosity and H 2 mass (aco = 3.5 M e (K km s ' 1 pc~ 2 ) -1 ) and 
restricting the HERACLES sample to the COLD GASS stellar 
mass range, the mean molecular Td ep i for HERACLES is ~1 Gyr, 
which is consistent with the COLD GASS estimate. 

As mentioned in Sect. 14.21 Td ep i of NGC 4569 also shows a 
hint of radial trend-reaching values of ~0.7-0.9 Gyr within a ra¬ 
dius of 0.5 kpc. These low values for the molecular Td ep i are also 
in line with the recent results bv lFaesi et al.l (12014l) for NGC 300 
on 250 p c sc ales, so slightly larger than but comparable to ours, 
and bv lLada et al.l (1201 Oh for the Mil ky Way GMCs, wh o found 
depletion times of -0.2 Gyr. After iFaesi et ail (1201 4l) studied 
the relation between molecular gas and SF in 76 Hu regions of 
NGC 300, they concluded that the short Td ep i arises because their 
analysis accounts for only the gas and stars within the youngest 
star-forming regions. These depletion times correspond to the 
timescale for SF to consume the gas reservoir in the star clus¬ 
ters’ parent GMCs, which may be the more relevant quantity in 
the context of GMC-regulated SF in galaxies. 

The hint of a radial trend of Td ep i in NGC 4569 also sug¬ 
gests that SF becomes more inefficient at larger radii under the 
assumption of a fixed Aco conversion factor. This finding agrees 
with results obtained bv lLerov et al. ( 2013) who find systemati¬ 
cally lower Tdepi for a given Aco in the inner kpc of a sample of 
30 galaxies, among AGN and starbursts. Since Aco is typically 
smaller when closer to galaxy nuclei (iSandstrom et al.l 1 20131) . 
implying shorter Td ep i for a given CO luminosity, the radial trend 
in Tdepi could be even steeper. The shortening of Td ep i toward the 
very center of NGC 4569 coincides both with a better correla¬ 
tion between In, and Isfr distributions (see Fig.Q} and with an 
increase in the /G 1 ratio there, which is indicative of more ex¬ 
cited gas. The value of R 21 ranges from 0.2 to 0.9 in the inner 
region of NGC 4569 dBoone et al .1120071 1201 ll) . with a smooth 
and fairly symmetric distribution with respect to the center, and 
continuously increasing toward the nucleus. 

Similar radial trends in /Githat are also present in other 
NUGA galaxies (NGC 4579, NGC 6574: iGarcia-Burillo et al.l 


Lindt-Rrie g et all 20081 respectively) and in the sam- 
Lerov et al.l ( 2013 ) presumably reflect similar changing 


12009 : 
pie of 

physical conditions that drive the lower Aco factors found by 
ISandstrom et all ( 1201 3b . This would underscore that molecular 
gas in the central regions gives off more CO emission, appears 
more excited, and forms stars more rapidly than molecular gas 
farther out in the disks of galaxies. Compared to /G 1 . R 32 shows 
a clumpier and mor e asym metric distribution within the central 
region dBoone et al.ll201 lb- differently from the large-scale gra¬ 
dient in R 32 reported bv lwilson et al.l (120091) . Usually, R 32 is ex¬ 
pected to probe more extreme physical conditions, such as those 
occurring in star-forming or shock regions. fBoone et all ( 2011 ) 
indeed confirm that R 32 tracks massive SF in NGC 4569 bet¬ 
ter than /Gi by comparing /Gi and R 3 2 maps t o the Paa- im¬ 
age as detected by NICMOS (see also IWilson eFal . 2009 ). In 
any case, in NGC 4569, 12 CO is less excited (R 32 ~ 0.6) than 
in the Galactic center and in most of the centers of normal, 
starburst, and active galaxies studied so far for which ~ 


0.2 - 5 (e.g., Devereux et al. 199 z 

1 Mauersberger et al. 

1999 

Matsushita et al. 2004; Mao et al. 

201 of Combes et al. 

2013 
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Table 6 : Stellar and star formation properties of our galaxy sample. 


Galaxy 

( 1 ) 

log(M.) (a) 
[M 0 ] 

( 2 ) 

log(TIR) (b) 

[ergs -1 ] 

(3) 

SFR ,c) 

[Mo y r_1 ] 
(4) 

log(sSFR) 

[yr -1 ] 

(5) 

NGC 3627 

10.57 + 0.13 

10.4 

4.33 

-9.93 

NGC 4569 

10.38 + 0.12 

9.7 

0.86 

-10.44 

NGC 4579 

9.96 ± 0.23 

10.1 

2.17 

-9.62 

NGC 4826 

9.99 + 0.12 

9.6 

0.69 

-9.30 


Notes. (a) Values from lSkibba et alJ l feOlllt . (i) Values from iDale et alj d2012l ) . <c) SFR derived adopting calibration from lKennicutti dl998ali for TIR 
luminosity. 


iGarcfa-Burillo et all 12014)) . Shorter Td ep i could also be due to 
the high pressure driven by the deep potential well in the cen¬ 
tral parts of galaxies, driving gas to higher densities. This effect 
has been seen both in our G alaxy (lOka et al.ll2001h and in others 
(iRosolowskv & Blitzl2005l) . and it has been suggested that it ex¬ 
plains shorter T^iJ^Le., enhanced SFE in galaxy centers) in the 
sample of lLerov et ali (120131) . 

NGC 4579 and NGC 4826 instead have higher molecular 
Tdepi of ~2 Gyr, both with I 2 CO(1-0) and 12 CO(2-l) when 
only taking the statistically significant cases into account. For 
NGC 4579, neglecting the case of the net sublinear K-S re¬ 
lationship found with 12 CO(2-l) at the resolution of 76 pc, 
Tdepi ~2 Gyr is accompanied by an index Vn t around unity. As 
already mentioned above, the interpretation of a linear K-S re¬ 
lationship is that CO is primarily tracing star-forming clouds 
with relatively uniform properties, including Esfr, and conse¬ 
quently the Tdepi of the CO traced gas is constant and of ~2 Gyr, 
both within and between galaxies. However, there is still no 
consensus on either the preci se K-S parameter estima tes or the 
associated interpretation (see iDohhs & Pringlell2013l for a re¬ 
view of explanations of the K-S relationship). In contrast to 
NGC 4569, Tdepi in NGC 4579 does not show a radial trend al¬ 
though R 2 i varies as a function of the radius. Close to the AGN, 
at r < 200 pc, Rj \ « 0.8-1, whereas at r > 500 pc /A i ~ 0.4-0 .6 
(IGarcfa-Burillo et al .1120091) . 

The use of the 12 CO(3-2) line gives a Td ep i of ~1 Gyr both 
for NGC 4569 and NGC 4826, suggesting that the gas traced 
by this transition could be more efficiently converted into stars 
with respect to 12 CO(1-0) and 12 CO(2-l) emissions, despite us¬ 
ing mean CO(l-0)/CO(3-2) ratios to convert mass densities. 
While NGC 4569 has Td ep i ~ 1 Gyr for all three 12 CO lines 
and for all resolutions taken into account, NGC 4826 shows a 
lower 12 CO(3-2) Td ep i than those derived from 12 CO(1-0) and 
12 CO(2-l) (~2 Gyr). In addition, for NGC 4826, this trend is 
seen by using 12 CO(3-2), while Td ep i does not show a radial 
trend in the cases of 12 CO(1-0) and 12 CO(2-l). However, the 
12 CO(3-2) distribution is similar to that of the CO(l -O) and 
CO(2-l) lines at the same resolution dBoone et al.l201 lh . In con¬ 
trast to NGC 4569, R 32 does not appear to be clumpier than 
R 21 in NGC 4826, and both maps look relatively symmetric 
with respect to the dynamical center, with both ratios increasing 
steadily toward the center (iBoone et alil201 ll) . As in NGC 4569, 
high R 32 values (>0.3) occur within the Paa emission region 
as detected by NICMOS, while this is not true for /G 1 , which 
reaches high values (>0.7) away from the Paa emitting regions 
dBoone et al.ll201 lh . This difference, along with the variation in 
R 32 , could be responsible for the short Tdepi inferred by 12 CO(3- 
2) in NGC 4826. 


Compared to the 12 CO(1-0) and l 2 CO(2-l), the 12 CO(3-2) 


traces relatively warmer and denser gas. There is also growing 
evidence that the 12 CO(3-2) emission correlates more tightly 
with the SFR or SFE than does 12 CO(1-0) and 12 CO(2-l) (e.g 
Komugi et alJ l2007t Murao taefalj 120071 [Wilson et alJ [2009, 
2012 ). Similarly, Iono et al. ( 20091) show that the 12 CO(3-2) 


emission correlates nearly linearly with the FIR luminosity for 
a sample of local luminous IR galaxies and high-redshift sub¬ 
millimeter galaxies. Thus, it appears that the 12 CO(3-2) emis¬ 
sion primarily traces the molecular gas associated directly with 
SF, such as high-density gas that is forming stars or warm gas 
heated by SF, rather than the total molecular gas content of a 
galaxy usually traced by 12 CO(1-0) and 12 CO(2-l). Neverthe¬ 
less, this evidence is not able to explain the shorter Td ep i char¬ 
acterizing the 12 CO(3-2)-molecular K-S relation for NGC 4826 
with respect to 12 CO(1-0) and 12 CO(2-l). 

Values of Td ep i derived for a given galaxy and within compa¬ 
rable FoVs are indeed expected to be similar for all 12 CO lines, 
by construction, since we are not deriving gas traced by 12 CO(3- 
2 ) but rather the total gas, assuming a R 32 line ratio, in addition to 
the Xco conversion factor (see Eq.[T|). The bottom righthand pan¬ 
els of Figs. [7] and [8] show that Eh, data points derived from the 
12 CO(3-2) line emission have typically lower Eh 2 values with 


respect to CO(l-O) and l2 CO(2-l) derivation, but this effect is 
not due to the choice of the 12 CO line but simply to the availabil¬ 
ity of 12 CO(3-2) data only for NGC 4569 and NGC 4826. It is 


possible that our result is simply due to the variation in R 31 that 
could not be considered in the mean value used for the derivation 
of E h , in Eq.[4] 

The radial trend of Td ep i as a function of physical conditions 
of gas, typically expressed in terms of line ratios, is not clearly 
defined. While NGC 4569 exhibits a Td ep i that decreases toward 
the very center, together with an increase in /G 1 , NGC 4579 has 
constant Td ep i as a function of the radius but higher R 21 in the 
central galaxy region. This difference could be explained by the 
(dis)agreements between gas and SFR morphologies for the two 
galaxies. Both NGC 4569 and NGC 4579 are centrally peaked 
in Esfr without a counterpart in Eh 2 (see Fig. [lj. However, 
NGC 4569 shows CO emission in the very center, but NGC 4579 
has neither gaseous nuclear emission nor a radial trend in the gas 
morphology (although the TGi line ratio decreases as a function 
of radius). It appears therefore that the radial trend of Td ep i as a 
function of CO line ratios also has to be accompanied by sim¬ 
ilar radial trends in both gas and SFR morphologies in a given 
galaxy. 

It is well known that gas depletion times are different in ex¬ 
treme starburst galaxies and merging systems in the local Uni¬ 
verse, as well as in submillimeter galaxies at high redshifts 
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(e.g., Kennicuttlll998a: Bouc he et al.ll2007l : lRiechers et al .1120071: 


Bothwell et alJl2010h . In these galaxies, gas depletion times are 
significantly shorter than the normal galaxy population, with the 
star formation surface density being an order of magnitude larger 
at fixed gas surface density (iGenzel et al.ll2010t) . More recently, 
it has emerged that the molecular Td ep i varies even within the 
population of nearby ‘normal’ star-forming galaxies and that this 
variation has been quan tified as a functio n of a range of funda ¬ 
ment al parameters. fSaintong e et al.l (1201 ll) . iBoselli et al.1 (120141) . 
and iBothwell et al] ( 201 4) all find that depletion timescales vary 
monotonically from ~2 Gyr to ~0.1 Gyr across two orders of 
magnitude in stellar mass. This finding also persists when us¬ 
ing a c onstant conversion fac t or and different convers ion factor 
values dBothwell et al.l 120141) . ISaintonge et al.l (1201 ll) also find 
that Tdepi correlates inversely with the specific SFR (SFR per unit 
stellar masses, sSFR). Table [6] collects the main stellar and star 
formation properties of our galaxy sample. In this table. Col. (1) 
gives the galaxy name. Col. (2) the logarithm of the stellar mass 
(M,), Col. (3) the logarithm of the total IR (TIR) luminosity. 
Col. (4) the SFR derived adopting the conversion of iKennicutti 
d!998al) for TIR (from 8-1000 pm), and Col. (5) the logarithm 
of the sSFR. NUGA galaxies have stellar masses in the COLD 
GAS range (log(M„) from 10 to 11.5 M e ) and sSFR toward the 
high end of the COLD GAS sample (log(sSFR) from -11.8 to 
-9.5 yr _1 ), so are consistent both with our results for short de¬ 
pletion times and the correlation that lSaintonge et all (1201 ll) find 
between Td ep i and sSFR. Nevertheless, the NUGA galaxies with 
the longest Td ep i (although only a factor of 2) have the high- 
estsSFR, although the values are within the scatter found by 
ISaintonge et al.l (1201 ll) . 

The dependence of Td eD i on the spati al scale has also been 
explored by other groups. Schruba et al. ( 2010 !) find a median 
Tdepi ~1.1 Gyr for M 33 (using the conversion factor we use 
here) with no dependence on type of region targeted on kpc 
scales, but below these scales Td ep i is a strong function of both 
the adopted scale and the type of region that is targeted (see also 
Komugi et al. 2005, for Tdepi even dow n to 0.1 Gyr in the ce ntral 
few hundred parsecs of local spirals). fSchruha et al.l (l2010h de¬ 
rived very long Td ep i (up to 10 Gyr) using small apertures cen¬ 
tered on CO peaks, but very short Td ep i (0.3 Gyr) with small 
apertures targeted toward Her peaks. This may not be surpris¬ 
ing, given the disparity in the bright Her and CO distributions in 
M 33 (see Fig. 1 in lSchruba et al.ll2010l) . similar to the inconsis¬ 
tencies between X gas and Xsfr morphologies characterizing some 
NUGA galaxies, suc h as NGC 4569 and NGC 4579 (see Fig|TJ. 
ISchruba et al.l (l2010h find that the SF relation in M 33 observed 
on kpc scales breaks down for aperture sizes of ~300 pc. These 
aperture sizes are larger than the breakdown scale of 80-100 pc 
found bv lOnodera et al.l (1201 Oh for the same galaxy. 


5.3. The dispersion of the data 

At the relatively high resolutions studied here, the Pearson corre¬ 
lation coefficient r con is generally low, and there is large scatter 
with values as low as ~0.3. Angular resolutions of ~l-2" cor¬ 
respond to ~20-200 pc at the distances of the sample galaxies. 
An inspection of the maps shown in Fig. Q] reveals that CO in¬ 
tensity peaks are offset from the Ho peak, often by more than 
~l-2". In addition, there is the possibility that such small beams 
do not include both the star-forming region and its counterpart 
molecular cloud, which is generally several tens of parsecs in 
size. In general, star-forming regions do not coincide with the 
molecular clouds that give birth to them. This is typically seen 
in spiral arms of galaxies, where the optical and molecular spiral 


10 “ 


= 210“ 2 


10 - 



Fig. 8: Comparison between (i) the K-S relation plot derived 
for all sample galaxies and taking all the 12 CO emission lines 
at the common resolution of 200 pc into account, as shown 
in Fig. [7] and (ii ) results from IB igiel et afl (120081) . The model 
predictions from iKrumholz et al.l (120091) are also shown with 
solar metallicity and clumping factor c - 1, 2, and 5; see 
IKrumholz et al.l fl2009l) for more details. Also shown (as a blue 
short-dashed line) is_the trend for star-forming galaxies and 
mergers found bv lGenzel et al.l (120101) using a common conver¬ 
sion factor aco (= 3.2 M e pc 2 (K km s 1 ) 1 ). The horizontal 
axis has been adjusted to our aco (= 3.5 M e pc“ 2 (K km s -1 )' 1 ) 
and our non-inclusion of helium^ We have also plotted the 
merger/starburst trend by IGenzel et al.l ( 2010 ) (with aco = 
1 M 0 pc“ 2 (K km s' 1 ) -1 ) as shown by a (red) long-dashed line. 


structures are offset. If only the molecular clouds are observed 
and not the star-forming regions which are physically coupled to 
them, the gas density would be expected to show an excess over 
the expected correlation. 

For a given galaxy, for a given 12 CO transition, and with 
only the statistically significant cases, r corr is quite stable (see 
NGC 4569) or it tends to decrease with resolution, although this 
tendency is weak (see NGC 4579 in 12 CO(2-l) and NGC 4826 
in I2 CO(1-0)). The latter trend is consistent with the result of 
a large r scatter on smaller scales found bv IVutisalchavakul et al.l 
(120141) when studying the SF relation in the Galactic plane. They 
found that the rank correlation coefficients between log(Sn,) 
and log(LsFR) increase from ~0.4 to ~0.8 in the resolution 
range from 33" to 20'. lYutisalchavakul et akl (12014|) also com¬ 
pared their results with some resolved extragalactic SF relations, 
available on a scale of about 200 pc (for instance, M 33 from 
ISchruba et al.ll2010i) . finding the same trend of the dispersion as 
a function of resolution for extragalactic data and a smaller scat¬ 
ter in the depletion time distribution for the Galactic plane than 
for the extragalactic sources. 

5.4. The K-S relation at high densities of molecular gas 

The range of Xu, over which we studied the K-S relation is 
~20-120 M q pc 2 . These densities are higher than those of 
20-50 M q pc 2 studied by IBigiel etTiTJ ( 2008 ) and other au- 
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thors. The circumnuclear regions examined here are subject 
to physical conditions that are different from the disk. Be¬ 
sides denser gas tr acers, such as HCN and HCO + , observed in 
these regions (e.g., Krin s et ahll201 li [Combes et al.l l2013. 2014; 


iGarcf a-Burillo et al.ll2014l) . rotation curves in the central kpc are 

different from disks, following roughly rigid-body rotation com¬ 
pared to nearly flat in the disks (e.g.. ISofue et alJ 12003 ). This 
would lessen the effect of shear, which works on the clouds, 
and may enhance the efficiency of the SF arising from grav¬ 
itational collapse. These effects could conspire to produce a 
systematic change in the K-S relation in the central regions 
of galaxies. Figure [8] shows the comparison between our re¬ 
sults for the whole sample at 200 pc resolution (as_ also dis¬ 
played in Fig. |7| and the findings from lBigiel et al.|(| 2008). Also 


plotted in Fig. [8] are the predictions from the iKrumhol z et alJ 
( 2009j) models for solar metall icity and various clumping fac¬ 
tors (see IKrumholz et akll2009l for more details). For the com¬ 
pariso n we have converted the Xh 2 measurements in lBigiel et alJ 
( 200 8) to the conversion factor we use here, but have not ad¬ 
justed the SFRs since the Kroupa IMF ( Kroup a 120011) used by 
them is similar to ours (s ee ICalzetti et all 20071) . Figure [8] also 
shows two trends found bv lGenzel et al. ( 201 (it ): the (blue) short- 
dashed line gives the common trend for star-forming galaxies 
(SFGs) and mergers/starburst using a single conversion factor 
(aco = 3.2 M q pc~ 2 (K km s -1 )^ 1 ); the (red) long-dashed 
line corresponds to the merger/starburst trend only (with o'co = 
1 M 0 pc~ 2 (K km s -1 ) -1 ). Both regressions have been adjusted 
to our o'co (3.5 M e pc -2 (K km s I ) —I ) and to our not including 
helium. 

We find that at densities as high as ~ 10 3 M fi pc~ 2 there is 
reasonable agreement with the predictions of iKrumholz et all 
d2009l) . However, with a common conversion factor, the central 
regions of the NUGA galaxies have higher Xsfr for a given gas 
column than would be expected from the models. They fill in the 
region between the mergers or high-zstarburst systems and the 
more quiescent SFGs studied by iGenzel et al.l (120101 ). assuming 
that the former require a lower value of o'co- The apparently high 
Xh, w ould also be consistent with the results of Sandstro m et alJ 
d2013l) . who find that a co in centers of galaxies relative to their 
disks can be as much as ten times lower. We could thus be 
slightly overestimating Xh 2 in our sample, although it is diffi¬ 
cult to determine by exactly how much . To assuage the offse t 
between mergers or starbursts and SFGs. lKrumholz et al. ( 2012 ) 
have devised a “universal SF law” that is based on volumet¬ 
ric densities (see also lGenzel et aklfeOlOl) . They argue that such 
a formulation is particularly important for resolutions of 10- 
100 pc, similar to what is achieved by our observations. In addi¬ 
tion, the galaxi es studied here are roughly c onsi stent with model 
predic tions of[Rrumholz & VlcKe^ d2005h and Ikrumholz et all 
( 200 §. 


It is difficult to consider the K-S relation for normal galax¬ 
ies at even higher densities. This would require that we observe 
with even better spatial resolution (possible with ALMA), but 
the GMCs would then be highly resolved. Thus we may not 
be able to observe the molecular cloud and its counterpart star¬ 
forming region within one beam, for galaxies - such as those 
studied here - at distances of ~5-20 Mpc (see Sect. 15.31) . It is 
possible that the K-S relation is fundamentally a correlation that 
is seen only when the ISM is averaged over a considerable area. 
As already discussed in Sect. I5T1 our work shows that the K-S 
relation undergoes a change at a spatial resolution of ~80-100 pc 
(~GMC scale); i.e., its slope starts to decrease sharply with re¬ 
spect to higher scales. The breakdown of the K-S relation could 
therefore occur at ~80-100 pc, as found bv lOnodera et al.l ( 2010 ) 


for M 33, even if in NUGA dense nuclei a net breakdown is still 
not visible on scales as high as ~20 pc. 


6. Caveats 

Here we examine possible caveats or systematic uncertainties 
associated with our analysis and results presented in the previous 
sections. These can arise from physical causes connected to our 
simplified assumptions, from the methodology we used to derive 
the various quantities, and/or from the adopted fitting method of 
data. 


6.1. Derivation of I.] 


'Hi 


It is well known that the Wo conversion factor is a source 
of systematic uncertainty in the estimate of X H ,. By adopting 
a constant value for Wo- we do not account for its changes 
with galaxy pro perties, such as physical gas conditions and 


metallicity (e.g.. Malone 


Glover & Mac ] 


;., Maloney & B]ack 19881; Bolattoetaljj2008; 
Low! 1201 U Magrini et all 1201 It IBolatto et akl 


2013). Estimates of “typical” values of Wo in the Milky Way 


and other spiral galaxies range from ~1.5 x 10“” cm~ 2 (K km 
s -1 )' 1 to ~4 x 10 2 ° cm~ 2 (Kkm s -1 ) -1 . 

CO line ratios are typically used to constrain the excita¬ 
tion conditions (temperature, density, and column density) of the 
line-emitting gas. Some studies have shown that these line ra¬ 
tios (mainly Wi) te nd to radially vary, with higher values toward 
galax y nuclei (e.g., Lin dt-Krieg et al. 2008; Garcfa-Burillo et al.l 
12009 ). and this effect presumably reflects similar changing 
of the ph ysical condition s that drive the lower Wo factors 
found by ISandstrom et al.l (120131 ). When presenting a review 
of the theoretical underpinning, techniques, and results of ef¬ 
forts to e stim ate Wo conversion factor in different environments, 
IBolatto et al.l (i2013h also stress that Wo appears to drop in the 
central, bright regions of some but not all galaxies. 

The AGN may also play a role since a higher R 2 l toward the 
center of active galaxies could be interpreted as evidence of ex¬ 
ternal heating of molecular clouds by X-rays in the vicinity of 
the AGN (e.g JBaker et al.ll2003l ). However, this radial trend of 
the gas physical conditions does not unambiguously characterize 
all the LLAGN in our sample. Although R 22 ratios show a radial 
trend, or hints of it, in NGC 4569, NGC 4579, and NGC 4826, 
R 2l in NGC 3627 and R 32 in NGC 4569 and NGC 4579 do not 
exhibit the same behavior when assuming a constant or stochas¬ 
tic variation as a function of radius. For the lack of a com¬ 
mon trend in our sample, we adopted constant values of R 2 \ 
and R 22 . This assumption may introduce an uncertainty into the 
Xh, derivation, although NUGA galaxies benefit from individual 
measures of line ratio s for each galaxy nuc lear region. _ 

As for metallicitv. lBoissier et al.l(l2003l) and Naray anan et al.l 
(120 1 2l) find that by using metallicity-dependent Wo> £h 2 should 
typically change by a factor of 2 to 3 in the central regions com¬ 
pared to disk regions of nearby spirals, in the sense of decreas¬ 
ing the gas mass. Our sample galaxies have all been observed 
by Moustakas et al.l (2 010I) . who find roughly solar metallicities 
for all but were unable to evaluate metallicity gradients for lack 
of data. Nevertheless, because of the relatively small FoV af¬ 
forded by our CO observations, any metallicity gradient is ex¬ 
pected to be very slight over the radial range sampled by our 
data. We therefore adopted a constant Wo f° r the sake of unifor¬ 
mity and simplicity, as suggested in the review of iBolatto et al.l 
(120131) in the absence of a detailed characterization of Wo in a 
given galaxy and as typically done in the major portion of studies 
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aimed at determining molecular gas mass from 12 CO line inten¬ 
sities. 

As already discussed in Sect.|T| we only used molecular (hy¬ 
drogen) gas in our study based on the assumption that SF for¬ 
mation occurs in molecular cloud cores where FF dominates. 
This implies a two-step process of SF where H 2 clouds are first 
formed from Hi, and stars are then formed from H 2 . This leads 
to the hypothesis that H 2 should be more directly connected 
with SFR. Moreover, we are interested in galactic centers that 
are high-d ensity regions of galaxies and t hat ca n be Hi defi¬ 
cient (e.g., Bigi el et all 1200811 . lHaan et all (120081) observed the 
atomic gas for the galaxies in our sample, but the Hi beam ap¬ 
proximately subtends our entire CO FoV, ranging from ~20" to 
~40". With this beam, NGC 4579 has an Hi central hole with 
a surrounding ring, and this galaxy also has the lowest (Eh,) 
in our sample (see Table [4J. For the other galaxies, Hi column 
densities E ffl ~ 10-13 M 0 pc~ 2 for NGC 4569 and NGC 4826, 
and ~ 100 M Q pc -2 for NGC 3627. In the first two, mean Eh 2 
are 20-40 times higher, so neglecting Hi will introduce a few 
percent systematic underestimate of gas surface density; nev¬ 
ertheless, the Hi beam is large so this could be an underes¬ 
timate. NGC 3627 has the highest Eh, of our sample (600- 
700 M 0 pc 2 ), but the dense Hi gas is still a relatively small 
fraction. In any case, as shown by the model predictions of 
iRrumholz et all(l2009tl in Fig.[8j Eh, in the NUGA galaxies even 
at 200 pc resolution are in the surface-density regime where H 2 
is expected to dominate. 

For most of CO data of our sample, it has been possible to 
combine interferometric data with single-dish ones, allowing the 
“missing flux” problem to be resolved (see Table [3. Because 
in a interferometer two antennas cannot be placed closer than 
some minimum distance (Z) ni i n ), signals that are stronger than 
some value ( ocA/D m i „) will be attenuated. For only interferomet¬ 
ric 12 CO maps can the spatially extended sources therefore not 
be detected, which would have the effect of underestimating the 
gas surface density. 

6.2. Derivation of E S fr- diffuse Her, [N 11 ] contamination, and 
Ha extinction 

Our values of Esfr inferred from Hu could be overestimated if 
some fraction of the emission not directly ass ociat ed with ion¬ 
ization by massive stars were diffuse. iRahman et ah : (12011!) find 
that K-S slopes are steeper, the stronger the fraction of diffuse 
emission in the SFR tracer (see also lLiu et al.11201 1). However, 
they also conclude that diffuse emission either in the SFR tracer 
or in CO do not compromise results in the high column-density 
regions of the galaxy they studied, NGC 4254. Since NGC 4254 
has peak Eh 2 S200 M Q pc “ 2 , comparable to the mean Eh 2 of our 
sample galaxies, we do not expect the potential overestimate of 
SFR from diffuse Her emission to be a significant problem. 

We removed the [N 11 ] contamination within the filter band¬ 
pass using average [N n]/H a values available in the literature, 
as typ ically done in narrow-band studies. However. iBlanc et al.l 
( 2009 ) used integral-field spectroscopy to study the spatially- 
resolved SF relation in NGC 5194, and checked the validity of 
the assumption of a constant [N 11 \jHcr and the bias introduced by 
it. They find an increasing [N ii]/Hq' ratio as one goes to fainter 
Her fluxes, which is consistent with the nebular emission in the 
faintest parts of the galaxy (mainly in the interarm regions) be¬ 
ing dominated by the diffuse ionized gas. Assuming a constant 
[N n]/H a may thus introduce an underestimation of the Her flux 
(and therefore of the derived SFR) in the galaxy’s central re¬ 
gions. 


In SectionsfTland l3.3l we stressed that the Her emission lines 
are subject to extinction within the galaxies, which we took into 
account for our analysis by adopting Ay(int) values determined 
for each sample galaxy rather than a mean value. However, 
we ap plied this throughout, as suggested by iKennicutt & Kent] 
(119831) . This is somewhat of an oversimplification because cen¬ 
tral regions of galaxies are generally subject to greater extinction 
than disk regions. This may cause the extinction correction and 
(therefore) the SFR to be underestimated. 

6.3. Choice of the fitting method 

A central component of observational investigations is also the 
quantification of the correlation between two or more observed 
quantities. Typically, linear regression provides estimates of the 
zero point and slope of the “best-fit” regression line between 
the observed data. Since each adopted methodology for fitting 
a power law to the data introduces a bias in the results, it is im¬ 
portant to understand the limitations of common fitting meth¬ 
ods. Usually in the K-S relation studies a linear regression in 
logarithmic space is performed, but met hods differ in the treat¬ 
ment of error bars . Some works (e.g., IKennicutt et alJ 120071 : 
IRahman etak||20lU IMom ose et al.ll201 3l) used the FITEXY al¬ 
gorithm dPress et aDl 19891) . which has the advantage of incorpo¬ 
rating errors in both the ordinate and abscissa coordinates, thus 
providing robust regression results, although errors must be as¬ 
sumed to be symmetric in logarithmic space, which is not always 
the case. Other studies, already quoted in Sect. 13.71 and includ¬ 
ing the present work, adopted the OLS bisector method, which 
returns a bisector line in x- and y-axes without the errors being 
taken into account. Both methods have the disadvantage of not 
being able to incorporate upper limits in the minimization. 

For this reason, iBlanc et ah ( 200 9) introduced and used a 
new method for fitting, based on Monte Carlo approach com¬ 
bined with two-dimensional distribution comparison techniques, 
which is not affected by the above issue. This method has the 
strong advantage of including the regions not detected in the 
CO map, including those with negative measured fluxes, in ad¬ 
dition to incorporating the intrinsic scatter in the SF relation as a 
free parameter and performing the fitting in linear space, thereby 
avoiding the assumption of log-normal symmetric errors. For a 
detailed discussion of the disadvantages of the linear regressions 
with respect to a Monte Carlo fitting method for a more real¬ 
istic treatment of systematics and accurate determination of the 
para meters and errors, we refer the reader to IBlanc et al.l ([2009) 
and[Shettv et al] (120131) . Such an analysis is beyond the scope of 
this paper. 

Based on the IRAM PdBI CO data reduction performed for 
all NUGA galaxies, we preferred to take only data above 30- 
level into account, for which we are certain of the interfero¬ 
metric significance. When analyzing the SF relation at sub-kpc 
scales of ten nearby galaxies. IMomose et al.l(l2013l) note that the 
choice of CO threshold surface density (3cr also in their case) 
affects the results very little; however . when repeating the anal¬ 
ysis with various thresholds, they obtained the same findings in 
terms of the K-S relation. Since we had compared our results 
with lBigiel et al.l (120081) (see Fig. [8j who used the OLS bisector 
method, we also adopted it. 

We note that results derived from the robust regression-fitting 
method can differ very much from OLS bisector ones (see Ta- 
ble0]i, c onfirming that the choice of the fitting method affects the 
findings (IBlanc et al.l2009tlShettv et al.l2013l) . This further rein¬ 
forces the need to compare K-S relation results obtained based 
on the same fitting method. 
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7. Conclusions 


We have investigated the resolved molecular K-S relation in the 
central kpc of four NUGA galaxies that have different types, 
which interact or do or not, which are barred or not, which have 
either a Seyfert or a Liner nucleus. Their SFR ranges from -0.7 
to -4 M 0 yr~ 1 , and sSFR from 0.04 to 0.5 Gyr 1 . The spatial 
scales studied in this work lie between —20 and 200 pc, and den¬ 
sities sampled are higher than in the usual galaxy disks, extend¬ 
ing toward 100-1000 M 0 pc~ 2 . 

The first result was that each galaxy has its own molecular 
SF relation on the investigated spatial scales. The K-S relation 
is most often underlinear (see NGC 4569), with slopes ranging 
from -0.5 to -1.3. The derived depletion time scales range be¬ 
tween 1 and 2 Gyr, which is very compatible with what is found 
on a larger scale in the COLD GASS sample of nearby galax¬ 
ies. These results are valid regardless of the 12 CO line observed, 
from the (1-0) to (3-2) transition. We also found that the K-S re¬ 
lation changes behavior on spatial scales of -80-100 pc, with the 
index Wit starting to sharply decrease with finer resolution with 
respect to larger scales. A break in the relation between Eh, and 
Esfr is expected on small scales (-10-50 pc), owing to phenom¬ 
ena such as star formation feedback, the life time of clouds, a tur¬ 
bulent cascade, or magnetic fields ( e.g.. iRruiissen et akll2014l) . 
and it has_been observed in M 51 dLiu et al.l 120 lTF> and M 33 
(lOnodera et alJl20loh in dense nuclei. However, in dense nuclei 
available thanks to NUGA sample a net break is still not vis¬ 
ible on scales of 20 pc. This might be due to the higher den¬ 
sity of the GMCs present in galaxy centers, which have to re¬ 
sist higher shear forces. Finally, one of the most important find¬ 
ings was that the proportionality between Eh 2 and Esfr found 
between 10 and 100 M 0 pc~ 2 continues at densities as high as 
-10 3 M e pc~ 2 . However, with a common conversion factor, the 
central regions of the NUGA galaxies have higher Esfr for a 
given gas column than those expected from the models, with be¬ 
havior between the mergers or high-’ starburst systems and the 
more quiescent SFGs, assuming that the former require a lower 
value of aco- 
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